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Abstract 

Heat  transfer  experiments  were  conducted  using  a  heat  exchanger  behind  a  pulse 
detonation  combustor  and  a  Garrett  automotive  turbocharger  at  the  Air  Force  Research 
Lab  (AFRL).  The  equivalence  ratio  and  purge  fraction  were  held  at  1.0  and  0.9, 
respectively,  while  the  frequency  of  operation  was  varied  from  10  to  12  Hz  in  1  Hz 
movements,  and  the  fill  fraction  was  varied  from  0.5  to  0.8  in  0.1  increments. 
Temperatures  were  calculated  using  an  energy  balance  and  used  to  determine  turbine  exit 
enthalpy.  The  representative  turbine  inlet  enthalpy  was  calculated  using  compressor 
work  and  radiation  from  the  turbine.  Turbine  inlet  and  exit  temperatures  were  also 
measured  directly  using  J-type  and  K-type  thermocouples  and  compared  to  calculated 
values  using  the  heat  exchanger  approach.  Compressor  and  turbine  work  was  computed 
and  compared  with  recently  attained  values.  Efficiency  was  presented  for  varying 
pressure  ratios.  The  efficiency  measurements  were  compared  with  time  accurate 
efficiency  measurements  from  on-going  work. 
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ENTHALPY  MEASUREMENTS  OF  A  PULSE  DETONATION  DRIVEN 


COMBUSTOR 


I.  Introduction 


1. 1  Motivation 

In  recent  years,  Pulse  Detonation  Engines  have  gained  attention  for  the  promise  of 
improved  performance  over  conventional  turbine  engines  and  ramjets  (Dyer  2002).  More 
recently,  government  agencies  have  pushed  for  a  Pressure  Gain  Combustion  (PGC) 
engine  to  be  integrated  into  a  production  worthy  aircraft.  The  end  goal  of  this  program  is 
to  design,  build  and  test  an  engine  capable  of  powering  an  aircraft  through  various  flight 
regimes.  Due  to  their  promise  to  provide  increased  specific  thrust  and  decreased  fuel 
consumption  at  higher  speeds  than  conventional  Gas  Turbine  Engines  (GTEs),  PDEs 
prove  to  be  a  promising  portion  of  the  solution. 

The  solution  will  likely  take  the  form  of  a  hybrid  engine  that  will  incorporate  the 
standard  elements  of  a  GTE  (compressor,  turbine  and  nozzle)  and  replace  the  combustion 
chamber,  where  deflagration  of  fuel  would  normally  occur,  and  replace  it  with  a 
detonation  chamber.  There  are  a  number  of  steps  that  need  to  take  place  before  such  a 
hybrid  is  successfully  developed. 

PDEs  obtain  their  increased  efficiency  by  means  of  detonation,  a  pressure  gain, 
near  constant  volume  combustion  process.  Conventional  gas  turbine  engines  bum,  or 
deflagrate,  fuel  through  approximately  constant  pressure  combustion  in  the  Brayton 
cycle.  A  PDE  utilizes  detonations,  which  offer  much  higher  pressures  at  the  site  of  fuel 
ignition,  generating  less  entropy  in  the  process  and  ultimately  translating  into  more 
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energy  being  extracted  from  the  fuel-air  mixture.  The  PDE  cycle  is  often  associated  with 
the  Humphrey  cycle.  The  Brayton  and  Humphrey  cycle  are  compared  in  Figure  1-1  later 
on. 

A  simplified  model  of  PDE  operation  breaks  the  process  into  three  main  phases: 
fill  fire  and  purge.  During  the  fill  phase  fuel  and  air  is  premixed  at  a  specific  equivalence 
ratio  and  injected  into  a  tube  which  is  immediately  sealed  at  one  end  upon  completion.  In 
the  fire  phase  an  ignition  source  initiates  deflagration,  or  burning,  of  the  fuel  at  the  closed 
end.  As  the  flame  passes  through  the  fuel  air  mixture  pressure  inside  the  tube  builds  until 
the  deflagration  to  detonation  transition  (DDT)  occurs.  At  this  point  a  coupled  shock 
wave  and  flame  front  travel  through  the  remainder  of  the  fiiel/air  mixture.  This 
detonation  wave  of  combustion  gasses  is  allowed  to  exit  the  tube  producing  thrust.  In  the 
purge  phase,  the  tube  is  filled  with  air  alone  to  provide  a  buffer  between  fire  phases  and 
to  aid  in  cooling.  This  process  is  repeated  in  a  cyclic  fashion,  often  at  high  frequency. 

At  AFRL  testing  has  shown  that  a  PDE  is  capable  of  powering  a  turbocharger  as  a 
means  to  self-aspirate  (Hoke,  et  al.  2002),  a  critical  step  in  the  process  of  producing  an 
airworthy  engine.  It  has  also  been  shown  that  a  Pulse  Detonation  Combustion  (PDC) 
powered  turbo-charger  will  produce  more  specific  power  than  a  Steady  Deflagration 
Combustor  (SDC)  given  similar  operating  conditions  (K.  P.  Rouser,  P.  I.  King  and  F.  R. 
Schauer,  et  al.,  Unsteady  Performance  of  a  Turbine  Driven  by  a  Pulse  Detonation  Engine 
2010).  However,  there  is  still  a  need  to  determine  the  efficiency  of  a  PDC  driven  turbine 
compared  to  a  SDC  or  GTE.  In  order  to  make  this  comparison  the  PDC  exhaust  flow 
must  be  further  characterized,  to  include  measurement  of  the  turbine  inlet  and  exit 
temperatures,  and  pressures. 
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1.2  Problem  Statement 


The  primary  parameters  that  are  necessary  to  determine  enthalpy  and  efficiency  across 
a  turbine  are  turbine  inlet  and  exit  temperatures,  pressures  and  mass  flow  rates.  Due  to 
the  extreme  and  unsteady  temperatures,  pressures  and  velocities  of  the  flow  in  a  PDC 
driven  turbocharger,  acquiring  experimental  measurements  of  the  combustion  products  as 
they  leave  a  radial  turbine  can  be  difficult.  Observed  detonation  wave  speeds  for 
hydrogen-air  (Schauer,  et  al.  2005)  average  1800  m/sec,  and  detonation  temperatures 
average  above  2000  K.  Purge  air  is  typically  subsonic  and  at  ambient  temperatures. 
Reliable  measurement  techniques  have  not  been  established  for  such  a  regime;  therefore, 
a  method  must  be  developed.  A  one-dimensional,  constant  volume  thermodynamic 
analysis  will  therefore  be  used  as  a  means  for  comparison  to  a  one-dimensional,  time 
accurate  method  being  developed  in  parallel  (K.  P.  Rouser,  P.  1.  King  and  F.  1.  Schauer, 
etal.  2011). 

1.3  Research  Objectives 

The  ultimate  goal  of  this  research  was  to  gain  an  understanding  of  the  efficiency  of 
a  PDC  driven  turbocharger  and  to  enable  comparison  against  other  (Ramjet,  Gas  Turbine 
Engine)  cycles.  This  research  specifically  focused  on  the  flow  downstream  of  a 
turbocharger  or  a  PDC  when  the  turbocharger  was  not  present. 

The  first  objective  was  to  build  a  heat  exchanger  capable  of  cooling  the  turbine 
exhaust  to  a  temperature  that  could  be  measured  accurately  on  a  time  averaged  basis 
using  thermocouples.  Currently  it  is  not  possible  to  measure  the  temperature  of  PDC 
exhaust  flow  while  operating  at  high  frequencies  and  mass  flow  rates  directly  using  a 
thermocouple  or  similar  probe-like  device.  Other  non-intrusive  temperature 


3 


measurement  teehniques,  like  optical  pyrometry,  have  been  applied  with  some  success, 
but  are  limited  to  higher  temperatures.  These  limitations  do  not  allow  for  the  flow 
temperature  to  be  captured  over  an  entire  PDC  cycle.  The  exit  temperature  of  the  heat 
exchanger  exhaust  gas  along  with  the  energy  removed  from  the  gas  by  the  heat  exchanger 
were  added  together  to  provide  an  average  temperature  of  the  flow  over  a  complete  cycle. 

The  second  objective  is  to  determine  the  pressure  at  the  inlet  and  exit  of  the  heat 
exchanger.  Ideally  two  state  variables,  temperature  and  pressure,  will  be  obtained,  which 
will  then  allow  enthalpy  to  be  calculated.  Sonntag  (Sonntag  1991)  defines  enthalpy  as 

H  =  17  +  PK  Eq.  2.2.1 

or  the  combination  of  the  internal  energy  with  the  pressure  multiplied  by  the  volume  of 
the  system.  For  an  ideal  gas,  PF  can  be  restated  as  RT,  or  the  ideal  gas  constant 
multiplied  by  the  temperature.  Therefore  for  an  ideal  gas,  enthalpy  becomes  solely  a 
function  of  temperature.  The  enthalpy  gives  the  measure  of  the  total  energy  of  the 
thermodynamic  system.  This  simplification  is  not  applicable  to  efficiency,  therefore  the 
pressure  at  the  inlet  and  exit  of  the  turbine  is  still  necessary.  The  third  objective  was  to 
determine  the  efficiency  of  the  PDC  driven  turbocharger. 


1.4  Methodology 

This  experiment  employed  an  application  of  the  first  law  of  thermodynamics, 
utilizing  a  quasi-steady  approach.  A  heat  exchanger,  situated  downstream  of  the  turbine 
reduces  the  temperature,  and  to  some  degree  the  pressure,  of  the  flow  to  the  point  where 
measurements  may  be  made.  The  assumption  was  that  the  flow  exiting  the  turbine  from 
both  the  detonation  and  purge  phases  were  able  to  mix  in  the  heat  exchanger,  and  the 
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mixed  outlet  temperature  was  measured.  The  mass  flow  rate  and  temperature  of  this 
colder  exhaust  flow  and  the  coolant  used  in  the  heat  exchanger  were  combined  to 
determine  the  inlet  energy  of  the  flow  and  consequently  the  temperature  of  the  PDC 
exhaust  turbine  gases. 

Starting  with  a  general  form  the  conservation  of  energy  equation: 


F  =  F 

^stored 


E  +E 

out  generated 


Eq.  2.3.1 


The  system  is  defined  as  the  heat  exchanger  itself  The  system  does  not  generate 
energy  and  it  is  allowed  to  achieve  steady  state  in  this  experiment,  so  the  net  change  in 
energy  stored  is  zero: 


=  E„. 


Eq.  2.3.2 


Furthermore,  for  a  counter-flow  heat  exchanger,  the  total  energy  change  into  and 
out  of  the  system  is  the  sum  of  the  energy  change  in  the  coolant  and  exhaust  gas: 


E  +  E  =  E  +  E 

^coolant  gas  ^coolant  gas 

in  in  out  out 


Eq.  2.3.3 


The  time  rate  of  change  of  energy  term  for  any  system: 
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-m 


y  out 


h+-Y^+gz 
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+q+W 


Eq.  2.3.4 


The  heat  exchanger  is  not  doing  any  work,  so  the  W  term  becomes  0.  The 
specific  enthalpy  is  a  combination  of  the  thermal  energy  and  flow  work.  The  flow  is  not 
moving  vertically  so  the  potential  term  is  removed.  The  average  velocity  at  the  entrance 
and  exit  of  the  heat  exchanger  are  assumed  to  be  the  same  for  the  coolant,  so  the  term 
also  is  removed.  At  steady  state  conditions,  the  rate  of  change  of  energy  is  zero.  For 
ideal  turbine  gases,  the  change  in  specific  enthalpy  can  be  approximated  by: 


Eq.  2.3.5 


The  energy  rate  removed  by  the  coolant  water,  q,  can  be  determined  by: 


q  =  AH  =  m(C 


T 

p,in  in, coolant 


C  T 

p,out  out, coolant  ^ 


Eq.  2.3.6 


In  Eq.  2.3.7,  Tin  was  determined  using  the  temperature  of  the  coolant  entering 
and  exiting  the  heat  exchanger,  and  the  temperature  of  the  exiting  combustion  gases. 

This  calculated  inlet  temperature  provides  an  indicator  of  the  energy  in  the  flow  at  the 
turbine  exit.  The  T-s  diagram  in  Figure  I-l  further  illustrates  this  approach.  The  process 
(1-2)  is  isentropic  compression,  (2-3H)  is  constant  volume  heat  addition  or  detonation, 
(3H-4H)  is  isentropic  expansion. 
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Figure  I-l:  T-s  diagram  for  Humphrey  cycle  and  Brayton  cycle 

For  the  purpose  of  this  experiment  process  (1-2)  takes  place  in  the  facility 
compressor,  (2-3H)  in  the  detonation  tube,  (3H-3.5H)  across  the  turbine  and  (3.5H-4H) 
across  the  heat  exchanger. 

This  study  focused  on  3H  where  the  combustion  products  have  exited  the  PDC 
tube  and  3.5H  where  the  combustion  products  have  exited  the  radial  turbine.  The  fuel/air 
mixture  will  not  have  returned  to  ambient  conditions  and  will  still  have  retained  residual 
energy.  The  heat  exchanger  exhausted  into  the  test  cell  so  pressure  at  4H  was 
approximately  equal  to  pressure  at  1  (ambient).  The  energy  transferred  from  the  exhaust 
gases  into  the  coolant  was  measured  in  the  process. 

The  temperature  of  the  coolant  was  measured  at  the  inlet  and  outlet  of  the  heat 
exchanger.  The  PDC  exhaust  temperature  was  measured  at  the  outlet  of  the  heat 
exchanger.  The  inlet  temperature  of  the  heat  exchanger  is  a  function  of  the  change  in 
energy  of  the  coolant  added  to  the  residual  energy  left  in  the  gas: 
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T  = 

in, gas 


(wC  at)  +mC 

V  P  /water  ^ 


'  T 

p, gas, out  gas, out 


^gas  p, gas, in 


Eq.  2.3.7 
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II.  Background  and  Literature  Review 

Research  on  extracting  energy  from  the  flow  of  a  PDC  has  been  ongoing  at  the 
Air  Force  Research  Labs  for  almost  ten  years.  In  2002  it  was  demonstrated  that  a 
turbocharger  is  a  viable  method  for  self-aspirating  a  PDC  (Hoke,  et  al.  2002).  During  this 
experiment  the  PDC  driven  turbocharger  was  able  to  run  for  25  minutes  without  any 
noticeable  performance  reduction  or  visible  signs  of  damage.  Further  experimentation 
showed  the  ability  of  a  turbocharger  to  aspirate  a  PDC  while  still  producing  thrust  under 
varying  subsonic  conditions  to  include;  varying  frequencies  (20,  30  and  40  Hz),  fill 
fractions  (1.0,  1.5  and  2.0),  compressor  flow  rates  (10-20  Ib/min)  and  compressor 
pressure  ratios  (1.05-1.73)  (Schauer,  Bradley  and  Hoke  2003).  In  this  work  a  series  of 
high  speed  pressure  transducers  were  used  to  determine  the  effect  of  the  addition  of  a 
turbine  on  the  detonation  and  blow  down  process.  It  was  shown  that  the  turbine  has  a 
damping  effect  on  the  shock  that  is  driven  by  a  detonation  wave. 

More  recent  work  conducted  at  the  Air  Force  Research  Lab  has  showed  that  when 
powered  by  a  PDC,  a  turbocharger  extracts  more  power  than  when  driven  by  constant 
pressure  combustion  (K.  P.  Rouser,  P.  I.  King  and  F.  R.  Schauer,  et  al..  Unsteady 
Performance  of  a  Turbine  Driven  by  a  Pulse  Detonation  Engine  2010).  Furthermore,  this 
work  also  showed  a  41.3%  improvement  in  specific  power  and  27.8%  improvement  in 
Break  Specific  Fuel  Consumption  (BSFC)  with  a  PDC  in  comparison  to  constant  pressure 
combustion.  It  is  important  to  note  that  this  increased  power  was  attained  at  low  pressure 
ratios.  This  magnitude  of  improvement  in  specific  power  is  not  expected  to  be  duplicated 
at  higher  pressure  ratios,  however  the  trend  of  increased  specific  power  from  a  PDC 
driven  turbocharger  versus  one  powered  by  SDC  is  still  anticipated.  To  reach  the 
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conclusions  on  improved  performance  high  speed  (5  MHz)  compressor  pressure,  mass 
flow  rate  and  tachometer  data  was  used.  Later  work  (K.  P.  Rouser,  P.  1.  King,  et  al., 
Parametric  Study  of  Unsteady  Turbine  Performance  Driven  by  a  Pulse  Detonation 
Combustor  2010)  at  the  Air  Force  Research  Lab  used  similar  measurement  devices  to 
show  that  the  average  specific  work  performed  by  a  PDC  driven  turbocharger  increased 
directly  as  a  result  of  higher  operating  frequencies.  At  those  higher  frequencies  rotor 
speed  response  approached  quasi-steady  behavior,  or  the  variation  between  the  peak  and 
minimum  rotor  speed  decreased  as  frequency  increased. 

Rouser  et  al.  also  evaluated  a  number  of  other  approaches  to  acquiring  flow  field 
data  of  a  PDC  driven  turbine  (K.  P.  Rouser,  P.  I.  King  and  F.  1.  Schauer,  et  al.  2011). 
More  specifically,  this  work  measured  turbine  rotor  speeds,  turbine  pressure  ratios  and 
flow  temperatures,  velocities  and  densities  of  the  unsteady  exhaust  leaving  a  PDC 
powered  turbocharger.  Future  applications  of  this  work  include  determining  unsteady 
turbine  efficiency.  Rouser  et  al.  also  noted  the  need  for  other  formulations  for  unsteady 
turbine  efficiency,  a  requirement  which  this  experiment  hopes  to  satisfy. 

PDC-turbine  integration  work  has  also  been  performed  incorporating  axial  flow 
turbines  and  an  array  of  PDC  tubes  with  a  bypass  ratio  of  7  (Glaser,  Caldwell  and 
Gutmark  2006),  vice  a  radial  turbine  powered  by  one  PDC  tube  as  in  previously 
discussed  works.  Recorded  turbine  inlet  temperatures  were  considerably  less  than  will  be 
reported  in  this  work,  mainly  due  to  the  quantity  of  bypass  air  that  was  used.  According 
to  this  study  increasing  fill  fraction  had  the  effect  of  increasing  turbine  inlet  temperature, 
specific  power  and  efficiency. 
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In  further  investigations  in  2007,  Glaser  et  al.  compared  the  efficiency  of  a  PDC 
driven  turbine  with  a  steady  flow  combustor  driven  turbine(Glaser,  Caldwell  and 
Gutmark  2007).  This  work  was  very  similar  to  that  of  Rouser  et  al.  with  the  exception 
that  the  radial  turbine  is  replaced  with  an  axial  turbine.  Bypass  flow  was  added  and  an 
array  of  six  tubes  was  used.  Glaser  et  al.  compared  the  turbine  efficiencies  of  the  two 
arrangements  by  using  turbine  inlet  temperature  and  pressure  ratio  across  the  turbine. 
Direct  measurements  were  possible  due  to  the  cooling  effect  of  the  bypass  flow.  Results 
showed  that  the  efficiency  of  the  PDC  driven  axial  turbine  was  comparable  to  the  steady 
combustor  driven  turbine. 

Rasheed  et  al.  also  experimented  with  a  single  stage  axial  turbine  using  60  % 
bypass  flow  in  2005(Rasheed,  Furman  and  Dean  2005).  This  work  reported  compressor 
work  of  100  hp  with  a  primary  mass  flow  rate  of  2  Ib/sec  (this  includes  fill  and  purge  for 
eight  tubes)  and  a  bypass  flow  rate  of  3  Ib/sec.  A  maximum  of  350  hp  was  calculated 
with  a  primary  mass  flow  rate  of  3  Ib/sec  and  a  bypass  flow  rate  of  5  Ib/sec  while 
operating  at  20  Hz. 

Further  analytical  work  has  included  resolving  the  efficiency  for  a  turbine  under 
unsteady  and  periodic  flow  conditions(Suresh,  Hofer  and  Tangirala  2009).  The  flow 
behind  a  PDC  would  match  this  description.  Suresh  et  al.  recognized  several 
formulations  for  this  efficiency,  one  variant  involving  no  averaging  and  another  involving 
averaging  of  an  equivalent  steady  flow.  In  both  formulas,  the  main  question  lies  in  the 
development  of  the  ideal  case  (the  denominator  in  the  efficiency  formula).  In  the  non- 
averaged  efficiency,  the  ideal  change  in  enthalpy  is  determined  by  assuming  the  mass 
flow  through  the  turbine  is  expanded  instantaneously.  The  second  variant  works  by 
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averaging  the  mass  flow  over  a  cycle  and  mass  averaging  the  total  temperature.  The  total 
pressure  is  then  work  averaged.  This  work  notes  that  the  latter  formulation  is  sensitive  to 
form  that  the  averaging  takes.  Using  CFD,  Suresh  et  al.  determined  that  the  work 
averaged  efficiency  produces  a  result  that  is  approximately  10%  higher  than  the 
efficiency  provided  by  the  non-averaged  equation. 
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III.  Test  Setup 


III.l  Preparation 

The  first  goal  in  this  experiment  was  to  determine  the  design  for  the  heat 
exchanger.  Utilizing  heat  transfer  principals  (Incropera,  et  al.  2007)  a  spreadsheet  was 
developed  to  iterate  on  heat  exchanger  design  parameters.  The  formulas  used  and  a  copy 
of  the  spreadsheet  can  be  found  in  Appendix  A.  This  spreadsheet  provided  the  initial 
estimates  for  heat  exchanger  dimensions  given  the  following  inputs: 


Table  III-l:  Heat  exchanger  spreadsheet  inputs  and  outputs 


Inputs 

Temp,  hot  gas  in 

1450  K 

2150  F 

Temp,  hot  gas  out 

900  K 

1160F 

Temp,  liquid  in 

293  K 

68  F 

Mass  flow,  hot  gas 

.04  Kg/s 

5.3  Ib/min 

Mass  flow,  liquid 

.35  Kg/s 

5.5  gal/min 

Inner  pipe  diameter 

.089  m 

3.5  inch 

Outer  pipe  diameter 

.162  m 

6.35  inch 

Outputs 

Temp,  liquid  out 

313.75  K 

105.35  F 

Heat  exchanger  length 

5.97  m 

19.6  ft 

Overall  heat  transfer 
coefficient 

26.7  W/m^K 

III.2  Data  Collection  Instrumentation 

Over  the  course  of  this  experiment,  the  test  setup  varied  significantly.  For  a 
detailed  description  of  the  previous  designs  that  led  to  the  final  experimental  setup,  see 
Appendix  D.  The  following  is  a  description  of  the  final  setup  that  was  used  to  produce 
the  results  discussed  in  Chapters  IV  and  V  this  paper. 
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A  counter  flow  heat  exchanger  design  was  chosen  due  to  its  increased  heat 
transfer  rate,  q,  per  unit  of  surfaee  area.  Schedule  10  (l/8th  ineh  thick)  6061-T6 
aluminum  was  selected  as  the  primary  building  material  for  its  favorable  eonductive 
properties  (180  W/m  K).  It  was  eventually  determined  during  the  course  of 
experimentation  that  one,  ten  foot  section  would  be  sufficient  to  reduee  the  temperature 
of  the  exhaust  produets  to  a  point  where  it  eould  be  directly  measured  with 
thermocouples.  The  three  inch  diameter  pipe  was  set  inside  of  a  six  inch  diameter  pipe 
and  a  ring  of  aluminum  was  welded  to  each  end  enelosing  an  annulus.  y4  inch  diameter 
holes  were  cut  in  opposite  ends  of  the  enclosure  facing  opposite  direetions  and  y4  ineh 
pipe  fittings  were  welded  to  each  hole  to  allow  for  coolant  flow. 


Figure  III-l:  Aluminum  heat  exchanger  section  before  addition  of  pipe  nipple 

Aluminum  pipe  nipples  were  welded  to  each  end  of  the  heat  exchanger  so  that  it 
could  be  seeured  to  the  exhaust  pipe  of  the  turbine.  Water  flowed  first  through  a  low- 
flow  liquid  flow  meter  upstream  of  the  heat  exchanger  shown  in  Figure  III-2: 
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Figure  III-2:  Low-flow  liquid  flow  meter 


The  coolant  then  flows  through  a  %  inch  diameter  Swagelok  fitting  and  across  a 
T-type  thermocouple  which  measures  the  inlet  temperature.  Next  the  water  flows  into  the 
annulus.  Flow  is  directed  vertically  into  and  out  of  the  heat  exchanger  in  order  to  reduce 
the  likelihood  that  air  pockets  form.  Water  flows  in  on  the  same  side  that  hot  gas  exits 
and  flows  out  where  the  exhaust  products  enter  the  heat  exchanger  (counter  flow  design), 
at  approximately  5.5  gallons  per  minute.  The  water  flows  out  through  a  similar  Swagelok 
fitting  where  another  T-type  thermocouple  measures  temperature.  Coolant  temperature 
measurements  are  taken  six  inches  prior  to  entering  and  six  inches  after  exiting  the  heat 
exchanger.  A  simplified  schematic  of  the  heat  exchanger  from  a  top  view  is  pictured  in 
Figure  III-3: 
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Figure  III-3:  Top  view  of  heat  exchanger 
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Pressure  was  measured  six  inches  upstream  and  downstream  of  the  heat 
exchanger  using  725  psi  static  pressure  transducers  with  a  one  foot  standoff  distance. 
This  standoff  distance  was  necessary  to  reduce  the  heating  of  the  pressure  transducer  and 
the  possibility  of  a  shock  wave  hitting  the  sensor  and  destroying  it.  Due  to  their 
sensitivity,  the  static  pressure  transducers  were  only  set  in  place  for  short  runs  (less  than 
ten  seconds)  to  further  reduce  the  risk  of  damage. 


Figure  III-4:  Sensotec  725  psia  pressure  transducer  fixed  to  detonation  tube  via  one 

foot  standoff  tube 

For  comparison  of  the  calculated  heat  exchanger  inlet  temperature,  during 
extended  runs  to  equilibrium  a  J-type  thermocouple  was  inserted  into  the  turbine  exhaust 
flow  at  the  location  where  the  pressure  transducer  was  during  the  shorter  runs.  J-type 
thermocouples  are  attached  on  the  exterior  of  the  heat  exchanger  three  inches  from  the 
inlet  and  the  exit  to  provide  temperature  measurements  used  in  calculation  of  the  amount 
of  heat  escaping  or  entering  the  heat  exchanger  via  radiation  and  or  natural  convection. 
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J-type  thermocouples  were  also  fixed  to  the  exterior  of  the  PDC  tube  at  25.5  inches  and 
48  inches  from  the  engine  head  for  the  same  purpose. 


Figure  III-5:  External  thermocouple  fixed  to  heat  exchanger 

An  additional  J-type  thermocouple  is  located  on  the  surface  of  the  aluminum 
nipple  used  to  attach  the  heat  exchanger  to  either  the  turbocharger  or  the  PDC  tube  in 
order  to  monitor  structural  integrity.  Figure  III-6  and  Figure  III-7  identify  the  locations 
of  thermocouples. 
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Figure  III-6:  Heat  exchanger  instrumentation,  head  end 


Figure  III-7:  Heat  exchanger  instrumentation,  tail  end 

The  experiment  was  earried  out  in  the  Pulsed  Detonation  Research  Facility  of  the 
Air  Force  Research  Laboratory,  using  a  similar  configuration  as  previous  work  (K.  P. 
Rouser,  P.  1.  King  and  F.  R.  Schauer,  et  ah,  LFnsteady  Performance  of  a  Turbine  Driven 
by  a  Pulse  Detonation  Engine  2010).  The  facility  supplies  compressed  air  to  the  main 
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and  purge  manifolds  as  seen  in  Figure  III-8.  The  PDRF  uses  Ingersoll-Rand  facility  air 
compressors  to  provide  compressed  air  to  the  manifolds.  Each  of  the  three  compressors  is 
capable  of  supplying  up  to  1412  ft  /min  air  mass  flow  at  pressures  up  to  100  psi. 
Compressed  air  flows  into  a  159  ft  receiver  tank,  and  then  is  routed  into  the  test  cell, 
where  it  is  separated  into  two  streams  for  the  main  and  purge  manifolds.  Each  air  stream 
is  controlled  by  Tescom  electromagnetic  controllers  that  actuate  pressure  regulators  and 
are  metered  through  calibrated  converging-diverging  nozzles.  Fuel  is  mixed  at  the 
entrance  to  the  main  manifold.  Fill  distribution  and  ignition  takes  place  using  an 
automotive  engine  head  and  cam  to  operate  intake  and  exhaust  valves  for  a  desired 
operating  frequency.  The  engine  head  is  taken  from  a  four-cylinder  engine  however,  for 
this  experiment  only  one  cylinder  head  was  used.  The  intake  valves  are  used  to  supply 
the  main  fdl  fuel-air  mixture,  and  the  exhaust  valves  are  used  to  inject  purge  air.  During 
the  fire  phase,  intake  and  exhaust  valves  are  closed.  The  turbocharger  and  heat 
exchanger  are  attached  downstream  of  the  engine  head,  or  to  the  right  in  Figure  III-8. 


Figure  III-8:  AFRL  Pulse  Detonation  Research  Facility  engine  test  block  diagram 


A  two-inch  diameter,  four  foot  long  steel  pipe  with  an  s-curve  was  used  for  a 
detonation  tube.  A  16  inch  Schelkin-like  spiral  assisted  the  DDT  process  by  increasing 


19 


the  turbulence  and  mixing  of  the  fuel  and  air  (Schauer,  et  al.  2005).  The  s-curve  was 
necessary  to  bring  the  exhaust  products  to  an  appropriate  height  to  make  Schlieren 
imaging  possible.  Three  ion  probes  are  installed  in  the  detonation  tube  to  verify 
Chapman-Jouget  velocities.  The  probes  short-circuit  when  the  flame  front  arrives,  and 
velocity  is  determined  from  the  transition  time  between  probes. 


Figure  III-9:  Two  inch  diameter  s-curve  detonation  tube  attached  to  engine  head 

Before  the  first  detonation,  the  turbocharger  turbine  is  driven  by  the  fill  and  purge 
phases  associated  with  the  start-up  sequence.  “Soft  starts”  were  used  for  all  runs  to 
prevent  detonating  a  larger  than  anticipated  volume  of  fuel  and  oxidizer.  In  the  soft  start 
process  the  ignition  source,  in  this  case  an  automotive  spark  plug,  is  initiated  prior  to  fuel 
being  added.  Fuel  is  gradually  added  to  the  main  air  until  the  desired  equivalence  ratio  is 
achieved.  PDC  operation  is  attained  by  first  setting  desired  operating  frequency, 
equivalence  ratio  and  fill  fraction.  Lab  View  software  determines  the  required  pressure 
to  achieve  the  given  fill  fraction  at  that  operating  frequency. 
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For  the  final  results  published  in  this  experiment  a  Garrett  GT2860RS  ball  bearing 
turboeharger  was  used.  This  turboeharger  uses  a  nine-blade,  radial  turbine.  The 
GT2860RS  is  also  equipped  with  a  radial  compressor  having  six  primary  impeller  blades 
and  six  splitter  blades.  The  GT2860RS  uses  a  76  trim  turbine  wheel  with  0.63  A/R 
turbine  housing  and  a  62  trim  compressor  wheel  with  a  0.60  A/R.  It  also  uses  a  T25 
turbine  inlet  flange  and  has  a  dual  ball  bearing,  oil  and  water  cooled  CHRA  (Center 
Housing  Rotating  Assembly). 


Compressor  Exit 


Turbine  Exit  Wastegate  Compressor  Inlet 


Figure  III-IO:  Garrett  GT-2860RS  ball  bearing  turbocharger 


The  turbine  inlet  of  the  turbocharger  is  coupled  to  the  PDC  exit  as  shown  in 
Figure  III-l  1  and  the  wastegate  is  disabled  so  that  all  of  the  mass  flow  from  the  PDC 
enters  the  turbine. 
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Figure  III-ll:  PDC  and  turbocharger  test  rig 


The  compressor  side  of  the  turbocharger  received  ambient  air  through  a  mass  air 
flow  (MAF)  sensor.  The  compressor  exhaust  pipe  consisted  of  a  two  inch  diameter  pipe 
and  a  ball  valve  to  back-pressure  the  compressor  with  a  50  psi  static  pressure  transducer 
and  J-type  thermocouple,  as  seen  in  Figure  III- 12.  Several  different  compressor  operating 
conditions  are  obtained  by  adjusting  the  PDC  operating  frequency.  J-type  thermocouples 
are  attached  to  the  turbine  housing  to  monitor  structural  integrity  of  the  turbine.  All 
temperature  data  is  sampled  once  per  second.  The  mass  air  flow  sensor  used  was  a  Pro- 
M  92mm  High  Flow  Mass  Air  Meter. 
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Figure  III-12:  Turbocharger  compressor  instrumentation  and  control  valve 


Figure  III- 13  illustrates  the  GT-2860RS  compressor  operating  map  where  target 
compressor  operating  conditions  run  down  the  center  of  the  efficiency  islands.  Operating 
frequencies  of  the  PDC  ranged  up  to  20  Hz  in  this  experiment,  with  fill  fractions  of  0.5 
through  0.8  and  purge  fraction  fixed  at  0.9.  The  fuel  being  detonated  was  hydrogen  with 
air  as  the  oxidizer  at  stoichiometric  conditions. 
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Figure  III-13:  Garrett  GT2860  compressor  operating  map 

The  exhaust  flow  of  the  heat  exchanger  needed  to  be  redirected  to  allow  for 
Schlieren  imaging.  To  accomplish  this,  a  six  inch  radius,  three  inch  diameter  180"  mild- 
steel  elbow  turned  the  flow  back  towards  the  engine  head.  To  orient  the  flow  in  the 
frame  of  the  Schlieren  camera,  an  additional  eight  feet  of  steel  pipe  was  used  for  the 
setup  without  the  turbocharger  and  nine  and  a  half  feet  was  used  for  the  configuration 
with  the  turbocharger. 


Figure  III-14:  180"  mild  steel  elbow  and  steel  pipe  extension  for  Schlieren  imaging 
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A  wooden  block  was  positioned  under  the  steel  extension  pipe  to  prevent  the 
elbow  from  yielding.  The  heat  exchanger  water  outlet  was  insulated  to  prevent  exhaust 
gases  from  affecting  the  measured  water  temperature. 


Figure  III-15:  Ten  foot  heat  exchanger  with  turbocharger 

Schlieren  imaging  took  place  at  the  exit  of  the  heat  exchanger,  with  the  purpose  of 
determining  the  density  gradients  in  the  exhaust  gases.  This  information  helped  to 
quantify  the  steadiness  of  the  flow  as  well  as  to  determine  relative  velocities  of  the 
exiting  flow.  The  Schlieren  system  in  Figure  III- 18  was  a  folded  z-type  arrangement 
consisting  of  12.5”  (3 1 .7  cm)  diameter  mirrors.  The  light  source  was  made  up  of  four 
pieces:  a  FSI  250W  halogen  illuminator,  two  50.8  mm  diameter  50  mm  focal  length 
lenses,  and  a  25  mm  by  25  mm  adjustable  slit.  The  camera  table  contained  a  vertical 
knife-edge: 
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Figure  III-16:  Phantom  camera  aimed  at  vertical  knife  (razor  blade)  edge 


Figure  III-17:  Phantom  ®  V7.1  high-speed  camera  with  a  l-2.8x  zoom  lens 


The  zoom  lens  allowed  focusing  onto  the  test  section. 
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Figure  III-18:  Schematic  of  Schlieren  setup 


The  Phantom  camera  could  record  a  maximum  of  75000  frames  per  second  at  256 
pixels  by  64  pixels  resolution.  The  camera  is  capable  of  higher  frame  rates,  and  higher 
resolutions,  but  not  both  at  the  same  time.  The  settings  used  are  a  compromise  between 
spatial  and  temporal  resolution. 

Schlieren  video  was  captured  at  10,  15  and  20  Hz  operation  with  and  without  the 
turbocharger  at  a  fill  fraction  of  0.5  and  a  purge  fraction  of  0.9  with  a  frame  rate  of  16000 
fps,  exposure  of  2  us,  resolution  of  256  by  256  pixels,  and  pixel  depth  (number  of  shades 
of  gray)  of  8  bit  (128  shades).  The  spatial  resolution  of  the  test  section  was  0.024  inches 
per  pixel  (.60  mm  per  pixel). 


III.3  Data  Collection  Procedure 

The  experiment  began  as  water  first  ran  through  the  heat  exchanger.  A  standard 
gate  valve  controlled  the  flow  rate  of  water.  The  water  passed  from  the  gate  valve  to  the 
low  flow  liquid  flow  meter.  This  flow  meter  used  a  12V  AC  power  source  and  outputs  a 
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sine  wave  with  a  K  factor  of  330  pulses  per  gallon.  The  flow  meter  output  is  connected 
to  an  oscilloscope.  The  frequency  (in  pulses  per  second)  is  recorded  and  then  multiplied 
by  the  K  factor  and  converted  to  minutes  to  determine  liquid  flow  rate  in  gallons  per 
minute. 
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=freq 
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Eq.  4.2.1 


The  gate  valve  was  manipulated  until  the  liquid  flow  rate  was  determined  to  be 
approximately  5.5  gallons  per  minute  via  Eq.  4.2.1.  The  actual  flow  rate  was  recorded  by 
hand. 

As  previously  mentioned  the  PDC  was  “soft  started”  during  all  experiments. 

During  the  soft  start  procedure  the  spark  was  started  and  fuel  is  gradually  added  to  the 
main  mixture  until  the  desired  equivalence  ratio  was  achieved.  The  detonations  are  then 
verified  using  the  ion  probes.  “Soft  starts”  were  necessary  to  prevent  fuel  and  air  from 
accumulating  in  the  heat  exchanger  and  turbine  and  then  being  detonated,  possibly 
damaging  both. 

One  of  the  key  test  parameters  that  were  varied  in  this  experiment  was  fill  fraction 
of  the  detonation  tube.  Fill  fraction  is  the  fraction  of  the  tube  that  is  filled  with  the  fiiel- 
air  mixture  prior  to  detonation.  Lab  View  software  in  the  test  facility  determines  the 
correct  pressure  and  mass  flow  rate  necessary  to  provide  the  demanded  fill  fraction.  It 
starts  by  determining  the  volume  of  tube  to  be  filled  by  multiplying  a  fill  fraction  by  the 
tube  volume,  both  of  which  are  input  by  the  user.  It  then  determines  mass  of  the  fill  by 
dividing  the  fill  volume  by  density.  The  fill  air  and  fuel  are  assumed  to  be  ideal  gases 
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and  density  is  considered  to  be  atmospheric  pressure  divided  by  the  gas  constant  and  the 
manifold  temperature.  The  mass  flow  rate  then  becomes  the  mass  multiplied  by  the 
frequency. 


Vr,„=(V...)(FF) 


The  software  calculates  mass  flow  rate  the  same  way,  regardless  of  obstructions  or 
manifold  pressures.  In  this  experiment  the  turbine  and  heat  exchanger  act  to  back 
pressure  the  system,  thus  increasing  the  amount  of  pressure  necessary  to  All  the  tube  to 
the  appropriate  volume.  Therefore,  without  further  adjustment  the  software  would  under 
All  the  tube  when  obstructions  were  present  downstream  of  the  engine  head.  For  this 
reason  the  ion  probes  are  oriented  so  that  they  straddle  the  cross  section  of  the  tube 
corresponding  to  a  All  fraction  of  0.5.  The  length  corresponding  to  a  All  fraction  of  0.5 
was  approximately  24  inches,  given  that  the  tube  was  48  inches  long.  The  estimated 
tube  volume  was  entered  and  the  tube  was  filled  to  what  was  thought  to  be  a  fill  fraction 
of  0.5.  A  short  run  was  performed  (just  long  enough  to  record  ion  probe  data).  If  the 
actual  fill  fraction  was  too  large  then  a  sharp  drop  was  noted  on  the  voltages  of  both  ion 
probes.  If  the  fill  fraction  was  too  small,  voltage  drops  were  not  observed.  To  correct 
either  issue,  the  tube  volume  parameter  in  the  Lab  View  software  is  manipulated  until  a 
sharp  drop  was  noted  on  the  first  ion  probe  and  no  voltage  drop  was  indicated  on  the 
second.  Using  this  technique,  it  is  possible  to  determine  fill  fraction  to  within  six  inches 
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(or  whatever  the  distance  is  between  ion  probes).  This  tube  volume  was  used  for  both 
configurations  (with  and  without  the  turbocharger)  in  order  to  ensure  that  a  similar  mass 
flow  rate  of  air  and  fuel  was  achieved.  This  fill  fraction  test  was  performed  before  the 
run  starts  to  ensure  the  appropriate  tube  volume  is  filled.  Fill  fraction  dissimilarities  are 
noted  in  section  IV.4  due  to  tube  heating,  which  caused  actual  fill  fraction  to  increase. 
The  tube  volume  used  for  this  experiment  was  168  cubic  inches. 

Two  separate  runs  were  used  to  gather  all  of  the  measurements  necessary  for  each 
set  point  (frequency  and  fill  fraction),  a  short  duration  run  to  collect  pressure  data  and 
Schlieren  imagery  and  a  long  duration  run  to  thermal  equilibrium  to  collect  temperature 
measurements.  For  the  short  duration  runs,  the  PDC  was  only  pulsing  for  approximately 
10  seconds,  just  enough  time  to  allow  for  pressure  measurements  to  be  captured  by  the 
transducers  and  for  the  camera  to  record  at  least  one  cycle  of  main  flow  and  purge  air 
leaving  the  heat  exchanger.  The  pressure  data  could  not  be  collected  continuously  up  to 
thermal  equilibrium  along  with  the  temperature  measurements  due  to  the  sensitivity  of 
the  pressure  transducers.  Ideally  the  pressures  would  have  been  recorded  at  thermal 
equilibrium  as  well  as  the  temperatures.  After  the  pressure  data  and  Schlieren  imagery 
were  collected  for  the  given  test  condition  the  pressure  transducers  were  removed  and  the 
port  was  fdled  with  a  bolt  and/or  a  thermocouple,  to  prevent  any  flow  from  escaping. 

The  high  speed  data  consisted  of  only  pressure  measurements  which  were  sampled 
at  5  MHz  while  the  rest  of  the  data,  to  include  mass  flow  rate  and  temperature,  was 
sampled  at  1  Hz.  The  high  speed  measurements  were  triggered  manually,  while  the  low 
speed  measurements  were  recorded  continuously  for  the  duration  of  the  run.  Each  run 
was  recorded  electronically  on  a  hard  drive  and  analyzed  at  a  later  time.  The  pressure 
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transducer  output  was  collected  via  a  Lab  View  program,  a  sample  of  whieh  can  be  seen 
here: 


Figure  III-19:  1  second  data  collection  at  10  Hz,  0.5  fill  fraction  and  0.9  purge 
fraction  of  pressure  downstream  of  turbocharger  (green)  and  downstream  of  heat 

exchanger  (blue) 


This  “high  speed”  data  was  then  converted  to  a  text  file.  During  the  conversion, 
sampling  rate  was  reduced  by  ten  times,  to  .5  MHz.  This  was  done  to  reduce  file  size  and 
to  reduce  the  load  placed  on  the  eomputer  in  manipulating  and  proeessing  data  with 
Microsoft  Excel.  The  one  half  second  long  pressure  fde  was  then  converted  from  voltage 
to  pressure  in  pounds  per  square  inch  absolute  via  the  formula: 


psia= 


f  Voltage^^^^^,^,  ^ 

^  Ipsia  ^ 

lOV  ^ 

gain  J 

to-139mvJ 

lVoltage^,^^i,^,i„J 

Eq.  4.2.2 
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When  the  pressure  data  downstream  of  the  turbine  and  the  heat  exchanger  were 
collected  with  the  725  psia  pressure  transducers,  it  was  necessary  to  put  a  gain  on  the 
output  so  the  variations  in  pressure  are  noticeable.  In  this  case  a  gain  of  40  dB  (voltage- 
ratio  of  100)  was  used,  which  corresponds  to  the  “gain”  term  in  Eq.  4.2.2.  The  second 
fraction  is  the  amount  of  voltage  generated  by  one  psi  and  the  last  fraction  takes  into 
account  that  the  voltage  used  to  power  the  transducer  was  higher  than  the  reference 
voltage  of  10  volts  and  the  pressure  transducer  was  not  regulated. 

The  low  speed  data  which  was  sampled  at  1  Hz  was  output  from  the  Lab  View 
software  directly  to  a  spreadsheet: 


UipCio&rd  ■-  1  Ortl  "3  Alinnmtn!  "  Mumtue  HiMi  C*ll«  tailjrtij 


Figure  III-20:  Low  speed  data  collection  spreadsheet 

For  extended  PDC  operation  to  thermal  equilibrium,  the  PDC  was  allowed  to  run 
until  the  water  temperature  of  the  heat  exchanger  and  the  exhaust  temperature  of  the  heat 
exchanger  remained  constant  (within  three  degrees)  for  approximately  30  seconds,  the 


PDC  set  point  was  then  changed.  Such  a  condition  was  considered  thermal  equilibrium. 
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Times  to  reach  thermal  equilibrium  varied  from  seven  to  ten  minutes  after  initial  start  up 
and  approximately  five  minutes  when  changing  between  set  points.  In  either  case,  the 
data  used  for  the  analysis  consisted  of  the  average  of  the  final  ten  seconds  of  each  run  to 
thermal  equilibrium. 
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Figure  III-21:  Data  collected  during  run  to  thermal  equilibrium  when  changing 

from  previous  set  point 


During  initial  experimentation  of  extended  runs  to  equilibrium,  the  heat  exehanger 
and  PDC  tube  were  allowed  to  eool  baek  down  to  the  water  inlet  temperature  and  room 
temperature  respeetively.  During  this  eool  down  period  the  PDC  was  still  running  as  it 
would  during  an  aetual  test,  however  fuel  was  no  longer  being  dispensed  and  the  spark 
was  diseontinued.  This  proeess  allowed  the  PDC  tubes  and  the  turboeharger  housing  to 


33 


cool  as  they  tend  to  reach  temperatures  in  excess  of  1200  °F  during  extended  runs.  After 
cooling  down  the  PDC  was  started  again  and  this  process  was  repeated  for  another  set 
point. 

Later  the  PDC  was  run  without  allowing  for  cooling  between  set  points.  The  water 
temperatures  and  heat  exchanger  exhaust  gas  rose  to  the  same  value  (within  3  degrees)  as 
when  cooling  was  allowed.  The  data  analyzed  in  this  experiment  was  collected  without 
allowing  for  the  system  to  completely  cool  in  between  runs,  except  where  it  was 
necessary  to  change  frequency.  When  changing  frequency  it  was  necessary  to  turn  of  the 
ignition  source  and  fuel  otherwise  the  PDC  would  backfire.  This  shutdown  interval  was 
minimized  as  much  as  possible  to  avoid  cooling  and  increasing  the  time  necessary  to 
reach  equilibrium. 

During  all  experimentation  the  test  cell  door  was  sealed  for  safety  and  the  status  of 
the  experiment  was  observed  via  remote  cameras  and  measurement  devices.  Table  III-2 
summarizes  the  operating  conditions  of  the  experiment  for  both  configurations: 


Table  III-2:  Experiment  operating  parameters 


Freq.  (Hz) 

10,  11,  12 

Fill  fraction 

0.5,  0.6,  0.7,  0.8 

Purge  fraction 

0.9 

(p  (eq.  ratio) 

1.0 

Ignition  delay 
(ms) 

3.0 

water  flow 
rate  (gal/min) 

5.6 

24  test  runs  were  necessary  to  collect  all  of  the  data  for  the  configuration  with  the 
turbocharger,  half  of  which  were  to  thermal  equilibrium  and  half  to  determine  pressures. 
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Another  24  test  runs  were  used  to  collect  the  data  for  the  configuration  without  the 
turbocharger. 


III.4  Data  Analysis 

The  experiment  was  run  with  a  ten  foot  long  heat  exchanger  to  reduce  the 
likelihood  for  condensation  inside  the  heat  exchanger.  When  water  was  not  condensing, 
the  gas  leaving  the  heat  exchanger  was  considered  to  be  the  combination  of  both  the  main 
and  purge  flows.  Similar  to  Eq.  2.3.7,  the  energy  conservation  equation  can  be  solved  for 
the  average  gas  inlet  temperature: 


(T), 


+q^+q,  1  -(rf,C,T) 


/cool 


main+purge 


Eq.  4.2.3 
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Eq.  4.2.4 


Note  that  the  radiation  (qrad)  and  free  convection  (q^)  from  the  heat  exchanger 
was  also  accounted  for.  Work  on  the  compressor  side  of  the  turbocharger  can  be  used  to 
calculate  turbine  work: 
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Eq.  4.2.5 


»t=H,-H,=m(C,T,-C,T,) 

The  temperature  change  across  the  compressor  actually  lags  behind  the  pressure  change. 
To  derive  the  work  performed  by  the  compressor  more  accurately,  the  isentropic  relations 
for  an  ideal  gas  are  used: 


T  =T 
^2  ^1 


vPiy 


r-i 


Eq.  4.2.6 


where  Ti  and  Pi  are  room  temperature  and  pressure  for  this  experiment.  Substituting 
Eq.  4.2.6  into  Eq.  4.2.5  produces  the  formula  used  to  calculate  compressor  work: 
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Eq.  4.2.7 


Over  the  course  of  this  experiment  it  was  found  that  the  turbine  housing  reaches 
temperatures  in  excess  of  800  K  depending  on  the  set  point  of  the  PDC.  Due  to  the 
significant  heating  of  the  turbocharger,  it  is  necessary  to  account  for  radiation.  To 
estimate  radiation,  the  turbine  was  approximated  to  be  a  six  inch  radius  sphere,  with  an 
emissivity  of  0.8.  The  following  equation  was  used  to  calculate  radiation: 

q,„,=^c7^(T/-T^^)  Eq.  4.2.8 
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This  analysis  was  performed  at  each  set  of  conditions  (frequency  and  fill 
fraction).  The  resulting  average  inlet  temperature  to  the  heat  exchanger  was  considered 
to  be  the  turbine  exit  temperature.  This  value  was  compared  to  the  measured  temperature 
at  the  exit  of  the  turbine.  The  turbine  inlet  temperature  was  determined  by  adding  the 
work  of  the  turbine  and  the  radiated  energy  from  the  turbine: 


'^C'pTj^b  +  ^comp+qrad 

T1  out 

mrb  = - TT^ -  or 

in  ^ 


^3 


Eq.  4.2.9 


The  compressor  work  and  mechanical  efficiency  are  used  to  approximate  turbine 

work: 


W 
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turb 


^mech 


Eq.  4.2.10 


The  efficiency  of  the  turbine  was  the  ultimate  target  of  this  experiment.  The 
efficiency  of  the  turbine  describes  how  well  the  turbine  performs  when  compared  with  an 
ideal  turbine.  An  ideal  turbine  is  one  that  does  work  reversibly  and  does  not  generate  any 
entropy,  or  energy  that  is  not  available  to  perform  useful  work.  The  actual  work 
performed  by  the  turbine  will  be  less  than  what  would  be  performed  by  an  ideal  turbine 
due  to  mechanical,  frictional  and  other  losses.  Formulations  for  efficiency  exist  for 
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turbine  efficiency  in  a  steady  flow  that  have  not  been  applied  to  unsteady  flows  such  as 
those  resulting  from  a  PDC.  The  following  formulation  uses  the  ratio  of  the  actual  work 
performed  by  the  turbine  to  the  ideal  work  performed  if  the  flow  was  expanded 
isentropically: 


H3-H,  _ 

H3-H4S  ^avg(C3,3T3-C^4j4s) 


Eq.  4.2.11 


where  T4S  is  the  temperature  at  the  exit  of  the  turbine  after  isentropic  expansion  and  Cp4s 
is  the  corresponding  specific  heat.  For  an  ideal  gas  with  constant  specific  heat,  T4S  may 
be  found  by  noticing  that  the  temperature  ratio  across  the  turbine  is  a  function  of  the 
isentropic  compression  ratio: 


vP4y 


This  equation  can  then  be  solved  for  T4S  producing: 


T  =T 

^4s  ^3 


y-l 


Eq.  4.2.12 


Eq.  4.2.13 


When  Eq.  4.2. 1 1  is  combined  with  Eq.  4.2. 13,  the  final  efficiency  equation  used  for 
calculations  in  this  experiment  results: 
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Eq.  4.2.14 
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Figure  III-22  helps  to  display  the  points  in  temperature  and  pressure  that  are  being 
pursued  in  this  experiment  on  a  temperature-entropy  diagram. 
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III.5  Uncertainty  Analysis 


An  uncertainty  analysis  is  necessary  to  determine  the  degree  of  accuracy  of  the  data 
provided  by  an  experiment.  The  total  uncertainty  is  the  square  root  of  the  sum  of  the  bias 
uncertainty  squared  and  the  precision  uncertainty  squared: 


Total 
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2  2 

W  +  p 


Eq.  4.2.15 


The  bias  uncertainty  can  be  thought  of  as  how  far  off  a  measurement  is  from  the 
actual  value.  Holman  presents  the  following  formula  to  determine  bias  uncertainty: 
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(Holman  1989)  Eq.  4.2.16 


where  wr  is  the  uncertainty  and  wi,  W2,  Wn  are  the  uncertainties  of  each  element.  Due  to 
the  complexity  of  the  formulas  used  in  this  analysis,  the  bias  uncertainty  will  be  broken 
down  into  parts  for  ease  of  calculation. 

Representative  values  were  used  to  calculate  the  bias  uncertainty  for  this 
experiment.  These  representative  values  are  similar  to  the  data  that  was  taken  in  the 
experiment,  and  have  a  slightly  higher  error  value  so  that  the  error  reported  is 
conservative. 

To  determine  wr  from  Eq.  4.2.16,  all  of  the  elemental  uncertainties  from  the 
system  must  be  addressed.  Detailed  calculations  of  the  uncertainty  values  can  be  found 
in  Appendix  F. 
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The  precision  error  can  be  thought  of  as  how  well  the  measurement  agrees  with 
itself  when  the  actual  value  being  measured  is  not  changing  or  the  randomness  of  the 
measurement.  For  measurements  that  were  recorded  electronically  with  multiple  points 
this  is  the  standard  deviation  of  the  samples  multiplied  by  the  inverse  of  the  t  distribution, 
The  following  table  summarizes  the  uncertainty  values  for  each  variable  that  was 
analyzed: 


41 


Table  III-3:  Uncertainty  values 


bias  uncertainty 
original  units  (percentage) 

precision  uncertainty 
original  units  (percentage) 

Kulite  725  psia  pressure 
transducer 

±  36  psia  (250) 

±0.2  psia  (1.4) 

Water  outlet  thermocouple 

+1.0  K  (5.0) 

±  1.4  K  (7.0) 

turbocharger  housing 
thermocouple 

±2.5  K  (.75) 

±7.8  K  (2.3) 

turbocharger  exhaust 
thermocouple 

+4.7  K  (.75) 

±1.4  K  (0.22) 

compressor  outlet 
thermocouple 

±2.2  K  (0.67) 

±0.1  K  (0.1) 

heat  exchanger  exit 
thermocouple 

±2.2  K  (0.4) 

±1.8  K  (0.7) 

water  mass  flow  rate 

±.06  gpm  (1.0) 

±.25  gpm  (4.5) 

mass  flow  rate,  gas  ( m  gas) 

±1.17E-4  kg/sec  (0.4) 

N/A* 

Reynolds  number  (Re) 

±113.6  (0.7) 

N/A* 

friction  factor  (f) 

±5.5E-5(0.2) 

N/A* 

Nusselt  number  (internal)  (Nu) 

±0.089  (0.2) 

N/A* 

convection  heat  transfer 
coefficient  (internal)  (h) 

±.429(1.6) 

N/A* 

heat  exchanger  purge  inlet 
temperature  (Tpuraein) 

±4.06  K  (0.7) 

N/A* 

Rayleigh  number  (Ra) 

±55300  (1.6) 

N/A* 

Nusselt  number  (external) 
(Nu) 

±0.332  (1.1) 

N/A* 

convection  heat  transfer 
coefficient  (free  convection) 
(h) 

±0.052  W/mK  (1.1) 

N/A* 

heat  transfer,  convection  (q^ee 

conv) 

±23.3  (4.6) 

N/A* 

heat  transfer,  radiation  (q^ad) 

±192  (1400) 

N/A* 

heat  transferred  to  water 

(Qwater) 

±25.5  W  (0.1) 

N/A* 

average  specific  heat  (Cp,avg) 

±12.2  J/kgK  (1.1) 

N/A* 

turbine  inlet  temperature 

(Tscorr) 

±15  K  (1.5) 

N/A* 

ideal  turbine  exit  temperature 
(T4s) 

±24  K  (1.9) 

N/A* 

enthalpy  gas,  out  (hgasout) 

±714  W  (1.0) 

N/A* 

heat  exchanger  inlet 
temperature  (Tj^am+purgejn) 

±49  K  (3.7) 

N/A* 

compressor  work  (Wcomp) 

±60  W  (4.7) 

N/A* 

turbine  work  (wt) 

±1408  W  (60) 

N/A* 

turbine  efficiency 

±  44  ("44^ 

N/A* 

*  Indicates  an  extremely  small  value  or  does  not  apply 
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IV.  Discussion  and  Results 


The  data  was  compiled  and  analyzed  as  described  in  Chapter  III  and  the  following 
trends  were  noted. 

IV.l  Heat  Exchanger  Inlet  Temperature 


Utilizing  Eq.  4.2.3  and  Eq.  4.2.9  the  turbine  inlet  and  exit  temperatures  were 
calculated  and  are  displayed  on  the  following  graphs: 


Temp  vs.  Fill  Fraction  10  Hz 


■  T4,calc 
♦  T4,meas 
AT3,calc 
XT3,  meas 


Fill  Fraction 


Figure  IV-1:  Measured  and  calculated  turbine  inlet  and  exit  temperature  at  10  Hz 
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Temp  vs.  Fill  Fraction  11  Hz 


■  T4,calc 
♦  T4,meas 
AT3,calc 
XT3,  meas 


Fill  Fraction 


Figure  IV-2:  Measured  and  calculated  turbine  inlet  and  exit  temperature  at  11  Hz 


Temp  vs.  Fill  Fraction  12  Hz 


■  T4,calc 
♦  T4,meas 
AT3,calc 
XT3,  meas 


Fill  Fraction 


Figure  IV-3:  Measured  and  calculated  turbine  inlet  and  exit  temperature  at  12  Hz 
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In  figures  IV- 1  through  IV-3,  the  T4,meas  was  the  average  temperature  recorded  by 
the  thermocouple  at  the  exit  of  the  turbocharger  and  T4  caic  was  the  temperature  calculated 
via  Eq.  4.2.3.  First,  it  should  be  noted  that  for  each  frequency  and  fill  fraction  the  turbine 
inlet  temperature  was  calculated  via  Eq.  4.2.9,  by  adding  the  work  performed  by  the 
compressor  and  the  heat  radiated  from  the  turbine  housing.  It  follows  that  the  changes 
from  T3  to  T4  measured  and  calculated  are  the  same  for  the  given  frequency  and  fill 
fraction. 

Due  to  the  lack  of  data  present  in  the  research  community  in  the  area  of  PDC 
driven  turbochargers  it  was  difficult  to  find  a  direct  comparison  to  the  values  derived  in 
this  research.  Most  non-intrusive  measurement  techniques  like  optical  pyrometry,  that 
capable  of  resolving  temperatures  at  high  speed  can  only  measure  temperatures  above  a 
certain  threshold.  For  this  reason  the  fire  phase  of  the  detonation  is  the  only  portion 
where  temperatures  were  recently  found.  Based  off  of  available  data  (K.  P.  Rouser,  P.  I. 
King  and  F.  I.  Schauer,  et  al.  2011)  the  average  temperature  at  the  turbine  exit  flow 
during  the  blowdown  phase  is  approximately  1600  K.  This  data  was  recovered  at  a  fill 
fraction  1.0  and  purge  fraction  of  0.5  at  15  Hz.  This  average  is  only  over  0.6 
milliseconds,  which  is  1/1 O**'  of  the  total  time  of  one  cycle  at  15  Hz  operation.  If  this  data 
is  averaged  with  a  conservative  estimate  of  temperature  during  the  rest  of  the  cycle  of 
1 000  K,  the  average  temperature  over  a  cycle  would  be  1 060  K.  This  temperature  falls 
very  closely  in  line  with  what  was  measured  via  the  thermocouple  at  the  exit  of  the 
turbine  during  the  most  demanding  operating  condition  (12  Hz,  0.8  fdl  fraction,  0.9  purge 
fraction).  Using  the  same  approach  for  the  conditions  at  the  inlet  of  the  turbine  an 
average  temperature  of  1 125  K  can  be  estimated.  This  temperature  is  also  very  close  to 
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the  turbine  inlet  temperature  calculated  based  off  of  the  measured  exit  temperature  of 
1092  K. 


A  number  of  trends  are  evident  in  figures  IV- 1  through  IV-3.  First,  there  is  an 
obvious  disparity  between  the  measured  and  calculated  inlet  temperatures  (T3,  meas  and 
T3,  calc).  The  measured  turbine  inlet  temperature  is  expected  to  be  higher  than  the 
calculated  turbine  inlet  temperature  due  to  minor  energy  losses  that  were  not  accounted 
for,  but  not  by  such  a  large  margin.  This  margin  is  most  significant  at  the  lowest  fill 
fraction  of  0.5.  This  margin  is  mostly  due  to  the  relative  speeds  of  the  different  gases 
and  the  fact  that  they  were  not  well  mixed  during  the  runs.  As  the  fill  fraction  increased 
the  temperature  calculated  at  the  inlet  of  the  heat  exchanger  agreed  more  with  the 
thermocouple  measurement  at  the  inlet.  This  is  due  to  the  fact  that  as  the  fill  fraction 
increased,  the  purge  fraction  remained  the  same.  At  a  larger  fdl  fraction  the  blow  down 
event  makes  up  a  larger  portion  of  the  exhaust  flow.  The  purge  fraction  at  this  larger  fill 
fraction  has  less  time  to  oscillate  and  is  in  contact  with  the  thermocouple  for  a  shorter 
period  of  time  in  comparison  to  the  same  purge  fraction  at  the  smaller  fdl  fraction.  This 
results  in  more  heat  from  the  blow  down  gas  being  transferred  to  the  thermocouple  and 
less  from  the  purge  gas  and  an  overall  more  accurate  temperature. 

Via  the  Schlieren  imagery  produced  in  this  experiment,  it  was  observed  that  the 
flow  is  still  segregated  at  the  exit  of  the  heat  exchanger.  The  turbine  helps  to  mix  the 
flow,  but  even  after  passing  through  the  turbine  ,  it  is  evident  that  the  hot  gas  is  moving 
much  faster  and  is  present  in  the  frame  for  a  much  shorter  period  of  time  (approximately 
1/6*  of  the  cycle  at  10  Hz).  The  Schlieren  videos  also  showed  that  cooler,  ambient  air 
was  actually  moving  back  into  the  heat  exchanger  in  between  fire  phases  due  to  the  low 
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pressures  in  the  tube  at  the  exit  of  the  heat  exehanger  eaused  by  the  blow  down  gas  as  it 
exited  the  heat  exehanger.  As  the  pressure  in  the  tube  equalizes,  the  purge  flow  expanded 
across  the  heat  exchanger  and  the  extension  pipe.  This  expansion  in  combination  with 
the  faster  moving  hot  gas  results  in  the  lower  temperature  purge  flow  contacting  the 
thermocouple  for  a  longer  period  of  time  in  comparison  to  the  blow  down  gas.  This 
causes  the  thermocouple  to  record  a  lower  overall  temperature  than  what  is  thought  to  be 
the  average  of  the  two.  Figure  IV-4  shows  the  density  gradient  at  the  exit  of  the  heat 
exchanger.  From  this  imagery  it  is  evident  that  the  flow  is  being  pulled  in  from  above 
and  below  the  pipe  as  is  indicated  by  the  density  gradients  curving  around  the  edges. 

From  watching  the  video  at  reduced  speed  it  also  becomes  apparent  that  the  hot  gas  from 
directly  in  front  of  the  pipe  is  also  being  pulled  back  in. 


Figure  IV-4:  Schlieren  images  indicating  suction  at  exit  of  heat  exchanger 
while  operating  at  10  Hz,  0.5  fill  fraction  and  0.9  purge  fraction 

Looking  at  the  differences  between  the  calculated  and  measured  temperatures  in 
figures  IV- 1  through  IV-3  it  becomes  apparent  that  the  higher  frequency  data  produces 
measurements  slightly  more  in  line  with  what  the  thermocouple  records.  The  increased 
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agreement  is  likely  due  to  the  faet  that  the  pressure  does  not  have  as  much  time  to 
oscillate  as  the  next  blow  down  phase  arrives  faster. 

If  the  pressure  oscillation  is  traveling  at  the  speed  of  sound,  which  is  approximately 
400  m/s  at  400  K,  it  would  take  it  approximately  0.015  seconds  to  reach  the  opposite  end 
of  the  heat  exchanger.  This  means  that  at  10  Hz,  with  no  obstructions  the  pressure 
oscillation  would  be  able  to  travel  up  and  down  the  heat  exchanger  approximately  five 
times  before  the  next  fire  phase  exits.  The  number  of  oscillations  is  actually  reduced,  as 
the  next  mass  of  exhaust  products  is  moving  towards  the  exit  and  the  pressure 
equilibration  process  taking  place  is  interrupted  by  the  next  blow  down  phase.  In  fact  at 
10  Hz  operation  and  a  fill  fraction  of  0.5  approximately  three  pressure  oscillations  are 
observed  in  the  Schlieren  video.  As  frequency  increases,  these  oscillations  have  less  of  a 
chance  to  reach  back  to  the  thermocouple  located  at  the  exit  of  the  heat  exchanger.  At 
approximately  65  Hz  operation,  the  pressure  pulse  would  not  be  able  to  travel  the  length 
of  the  steel  extension  pipe  back  to  the  thermocouple  in  time  to  beat  the  next  mass  of 
exhaust  products.  At  this  elevated  frequency  the  heat  exchanger  would  likely  become  a 
more  accurate  tool  for  measuring  inlet  temperature.  It  is  evident  from  the  results  above 
that  smaller  frequency  increases  help  slightly  with  accuracy. 

The  side  effect  of  reduced  suction  with  increasing  frequency  is  also  evident  in  the 
Schlieren  videos.  When  comparing  the  10  Hz  operation  to  the  15  and  20  Hz  operation,  it 
becomes  immediately  apparent  that  the  length  of  time  the  flow  is  retreating  into  the  tube 
decreases  as  frequency  increases.  When  comparing  the  exhaust  flow  of  the  PDE  driven 
turbo  and  the  PDE  alone,  it  is  evident  that  the  addition  of  the  turbocharger  reduces  the 
amount  of  time  the  flow  reverses  direction  and  moves  back  into  the  tube. 
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The  addition  of  the  turboeharger  also  has  the  effeet  of  mixing  the  flow.  Originally 
it  was  assumed  that  over  the  length  of  the  heat  exchanger  that  the  hot  gas  from  the  fire 
phase  and  the  cooler  flow  from  the  purge  phase  would  mix  significantly  inside  the  heat 
exchanger,  however  it  appears  that  without  the  turbocharger,  the  flow  is  still  considerably 
segregated.  The  stratification  of  the  flow  is  evidenced  by  the  darker  lines  in  the  Schlieren 
video,  signifying  larger  density  gradients  that  follow  immediately  behind  the  shock  as  it 
exits  the  exhaust  tube.  As  the  suction  occurs,  the  lines  become  lighter  as  the  tube  pulls 
cooler  air  back  inside.  The  thick  dark  lines  are  not  noted  again  until  the  next  shock  wave 
leaves  the  tube.  Suction  was  observed  at  the  exit  of  the  heat  exchanger  for  approximately 
.045  seconds  per  cycle  or  45%  of  the  cycle  time  at  10  Hz  operation  with  a  fill  fraction  of 
0.5. 


IV.2  Work  Calculations 


Compressor  work  was  calculated  from  Eq.  4.2.7.  All  values  in  this  calculation 
were  measured  directly  and  turbine  work  was  calculated  using  Eq.  4.2. 10.  The  change  in 
energy  across  the  turbine  was  also  calculated  as: 


r 


Eflux  ^turb 


A 


C  T  -C  T 

-^3  ^p  ^4 
corr  J 


Eq.  5.1.1 


where  T3  was  calculated  using  Eq.  4.2.9.  The  following  graphs  display  the  results: 
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Figure  IV-5:  Turbine  work  and  radiation  at  10  Hz  operation 
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Figure  IV-6:  Turbine  work  and  radiation  at  11  Hz  operation 
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Figure  IV-7:  Turbine  work  and  radiation  at  12  Hz  operation 


The  first  trend  noted  is  increasing  work  with  increasing  fill  fraction  and 
frequency.  This  trend  is  consistent  with  conclusions  from  other  work  at  AFRL  (K.  P. 
Rouser,  P.  1.  King  and  F.  R.  Schauer,  et  ah,  Unsteady  Performance  of  a  Turbine  Driven 
by  a  Pulse  Detonation  Engine  2010).  Based  off  of  this  work  it  would  be  reasonable  to 
achieve  increased  turbine  and  compressor  work  as  more  fuel  and  air  is  being  combusted 
at  higher  fill  fractions  and  frequencies.  Using  the  turbine  work  numbers  and  dividing  by 
mass  flow  provides  specific  work  which  is  a  better  unit  for  comparison: 
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Table  IV-1:  Specific  work  for  varying  fill  fractions  and  frequencies 


Fill 

Fraction 

Purge 

Fraction 

frequency 

10  Hz 

11  Hz 

12  Hz 

0.5 

0.9 

0.123 

0.153 

0.179 

0.6 

0.158 

0.179 

0.217 

0.7 

0.222 

0.246 

0.276 

0.8 

0.289 

0.307 

0.336 

1.0 

0.5 

0.288* 

0.334* 

0.368* 

*Reported  in  other  work 
R.  Schauer,  et  al.,  Unsteac 
Driven  by  a  Pulse  De 

K.  P.  Rouser,  P.  1.  King  and  F. 
Jy  Performance  of  a  Turbine 
Jtonation  Engine  2010) 

Comparing  the  values  of  specific  work  with  what  was  reported  recently  at  a  fill 
fraction  of  1 .0  and  purge  fraction  of  0.5  there  is  a  good  agreement.  Agreement  is 
expected  considered  the  values  were  measured  in  a  similar  fashion. 

One  might  expect  that  as  fill  fraction  and  frequency  increase  the  specific  work 
would  remain  relatively  constant  because  the  work  increases,  but  the  mass  flow  rate  is 
increasing  as  well.  This  increase  in  specific  work  is  likely  due  to  reduced  losses  at  higher 
rotor  speeds.  As  the  rotor  turns  faster  the  momentum  becomes  larger  in  proportion  to  the 
friction  acting  against  it. 

The  Eflux  term  is  substantially  larger  than  the  work  term,  and  this  is  mainly  due 
to  the  added  radiation.  It  would  be  expected  that  the  radiation  term  would  increase 
linearly  because  the  temperature  of  the  turbine  increases  in  a  fairly  linear  fashion  as  fill 
fraction  increases.  The  nonlinearity  of  the  turbine  Eflux  term  is  mainly  due  to  the  fact 
that  the  radiation  equation  utilizes  the  surface  temperature  raised  to  the  fourth  power. 
After  considering  this,  the  fact  that  the  Eflux  term  always  remains  approximately  50% 
above  the  turbine  work  term  makes  more  sense. 
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The  large  value  of  the  Eflux  term  would  indicate  that  the  turbine  cannot  be 
considered  as  adiabatic  or  the  efficiency  term  will  be  artificially  inflated  as  there  is  a 
large  portion  of  energy  that  is  being  leaving  the  flow,  but  isn’t  generating  any  work. 


IV.3  Efficiency 

The  overall  goal  of  this  experiment  was  to  help  determine  efficiency  for  a 
turbocharger  driven  by  unsteady  flow.  To  that  end,  both  the  measured  (via  the 
thermocouple  at  the  exit  of  the  turbocharger)  and  the  calculated  (via  the  heat  exchanger 
and  first  law  analysis)  turbine  exit  temperatures  and  their  respective  turbine  inlet 
temperatures  (calculated  based  off  of  turbine  work)  were  used  to  compute  a  notional 
efficiency.  To  compute  efficiency  the  pressure  ratio  across  the  turbocharger  was  varied 
and  the  rest  of  the  variables  were  held  constant  in  the  following  analysis.  The  actual 
pressure  data  that  was  recorded  using  the  725  psia  pressure  transducer  was  not  used  for 
this  analysis  due  to  the  magnitude  of  its  uncertainty  value.  All  efficiencies  were 
calculated  using  Eq.  4.2. 14.  For  each  fill  fraction,  the  same  turbine  work  term  was  used. 
The  charts  below  show  the  effect  of  turbine  pressure  ratio  on  thermal  efficiency  when 
turbine  temperature  ratio  and  mass  flow  rate  are  held  constant: 
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Varying  average  turbine  expansion  (P3/P4)  ratio 
vs.  efficiency  @  10  Hz  for  measured  turbine  exit 
temperature 
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Figure  IV-8:  Average  turbine  expansion  ratio  versus  efficiency  for  10  Hz 
operation  using  measured  turbine  exit  temperature  and  calculated  turbine  inlet 

temperature 


Varying  average  turbine  pressure  ratio  vs. 
efficiency  @  10  Hz  for  calculated  turbine  exit 
temperature 
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Figure  IV-9:  Average  turbine  expansion  ratio  versus  efficiency  for  10  Hz 
operation  using  calculated  turbine  exit  temperature  and  calculated  turbine  inlet 

temperature 
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Figure  IV-10:  Average  turbine  expansion  ratio  versus  efficiency  for  10  Hz 
operation  using  measured  turbine  exit  temperature  and  measured  turbine  inlet 

temperature 


In  Figure  IV-8  the  turbine  inlet  temperature  (T3)  used  to  calculate  efficiency 
resulted  from  Eq.  4.2.3  and  adding  the  turbine  work.  Figure  IV-9  the  efficiency  was 
calculated  using  the  turbine  inlet  temperature  from  the  measured  turbine  exit  temperature 
and  adding  the  turbine  work.  Finally  in  Figure  IV-10  the  turbine  inlet  temperature  used 
was  measured  at  the  exit  of  the  PDC  when  the  turbine  was  not  present. 

The  first  and  most  obvious  trend  in  the  efficiency  calculations  was  a  significant 
drop  in  efficiency  with  increasing  pressure  ratio.  At  first  this  seems  counter  intuitive  and 
that  increasing  pressure  ratio  should  increase  efficiency.  It  must  be  noted  that  this  is  for  a 
given  turbine  inlet  and  exit  temperature.  If  the  turbine  inlet  and  exit  temperature 
corresponded  to  each  pressure  ratio,  the  trend  in  the  charts  would  be  inverted,  that  is  to 
say  that  efficiency  would  increase  with  increasing  pressure  ratio. 

The  next  trend  observed  in  Figure  IV-9  is  that  as  fill  fraction  increases  from  0.5  to 
0.6  the  efficiency  decreases,  but  increases  from  0.6  to  0.8.  It  should  be  noted  that  this 
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variation  is  within  the  uncertainty  of  the  reported  efficiency.  The  overall  trend  noted  in 
the  three  figures  increasing  efficiency  with  fill  fraction.  This  trend  falls  in  line  with  the 
trend  of  increasing  specific  work  with  increasing  fill  fraction,  noted  in  the  previous 
section.  As  the  fill  fraction  increases,  the  turbine  will  see  an  increasing  mass  flow  and  as 
a  result  will  be  choked  for  a  larger  portion  of  the  cycle.  The  longer  the  turbine  is  choked, 
the  higher  the  efficiency  would  become,  until  it  is  choked  for  the  entirety  of  the  cycle 
where  efficiency  should  level  off. 

Recently,  time  accurate  data  was  recorded  for  pressure  across  a  PDC  driven 
turbocharger  at  15  Hz  with  fill  and  purge  fractions  of  1.0  and  0.5  respectively.  When  the 
inlet  and  exit  pressures  were  time  averaged  separately  and  then  divided,  the  result  was  an 
average  pressure  ratio  of  2.22.  Extrapolating  the  trend  of  increasing  fill  fraction  and 
applying  it  to  this  data  would  yield  an  approximate  efficiency  of  10%. 

For  a  T3  turbine,  the  maximum  efficiency  possible  is  approximately  72%  which 
occurs  when  the  flow  is  choked.  It  can  be  said  with  certainty  that  the  flow  through  the 
turbine  is  choked  at  a  pressure  ratio  greater  than  2.3.  The  time  accurate  data  is  presented 
below. 
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Figure  IV-11:  Time  accurate  turbine  inlet  and  exit  static  pressure  at  15  Hz 
with  fill  and  purge  fractions  of  1.0  and  0.5  respectively  over  one  cycle  (K.  P.  Rouser, 

P.  1.  King  and  F.  I.  Schauer,  et  al.  2011) 


The  average  pressure  ratio  of  the  time  accurate  data  indicates  that  more  than  50%  of  the 
time  the  turbine  is  unchoked,  meaning  the  efficiency  of  the  turbine  will  be  below  72%.  It 
is  also  apparent  from  the  time  accurate  data  that  for  a  portion  of  the  cycle  the  pressure 
ratio  is  actually  negative.  Taking  into  account  the  portion  of  time  the  turbine  is  unchoked 
and  negative  an  efficiency  of  10%  is  not  unreasonable. 


IV.4  Fill  Fraction  Disparities 


As  the  detonation  tube  temperature  increases,  the  density  of  the  fuel/air  mixture 
decreases.  The  decreasing  density  results  in  a  disparity  between  perceived  and  actual  fill 
fraction.  At  the  beginning  of  a  run,  a  voltage  drop  over  the  first  ion  probe  and  lack  of 
change  in  voltage  over  the  successive  ion  probes  were  noted  which  indicated  a  fill 
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fraction  of  approximately  0.5.  As  the  tube  temperatures  rose,  voltage  drops  across  the 
seeond  and  even  third  ion  probes  were  noted.  For  example,  a  test  run  began  and  ion 
probe  data  shows  the  fill  fraction  to  be  -0.5: 


Figure  IV-12:  Ion  traces  showing  fill  fraction  of  0.5  at  run  start 


As  the  detonation  tube  reaches  equilibrium  temperature  (approximately  680  F  for 
this  eonfiguration),  the  reactants  expand  at  an  increasing  rate  when  entering  the 
detonation  tube.  Therefore,  due  to  the  increased  temperature,  density  is  redueed  and  the 
reactants  occupy  an  increased  volume  of  tube  prior  to  detonation.  This  increased  volume 
registers  on  the  next  ion  probe  as  it  detonates,  implying  an  increased  fdl  fraetion: 
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Figure  IV-13:  Ion  traces  at  equilibrium  for  the  same  run  showing  increased  fill 

fraction 

This  phenomenon  is  unavoidable  for  extended  duration  runs  without  changing 
mass  flow  rate.  In  an  effort  to  reduce  the  variability  between  configurations,  mass  flow 
rate  was  kept  constant.  Due  to  the  fact  that  the  Lab  View  software  computes  fill  fraction 
based  off  of  the  tube  volume  input  by  the  user,  the  mass  flow  rate  of  the  main  and  purge 
flows  did  not  change  over  the  course  of  a  given  run.  Therefore,  the  user  need  only  to 
input  the  measured  tube  volume,  fill  fraction,  purge  fraction  and  frequency  and  the 
software  would  provide  the  necessary  mass  flow  rate,  regardless  of  back  pressure  placed 
on  the  system  (by  the  heat  exchanger  and/or  turbocharger). 

To  ensure  that  the  detonations  were  still  occurring  at  these  equilibrium  conditions, 
wave  speed  was  measured  at  equilibrium: 
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Figure  IV-14:  Ion  traces  at  equUibriuin  with  wave  speed  calculated  indicating 
Chapman-Jouget  velocity  for  hydrogen  and  air  at  stoichiometric  conditions 


For  hydrogen  and  air  at  stoichiometric  conditions,  typical  Chapman-Jouget 
velocities  are  approximately  1970  m/s(Schultz  2000).  For  the  drop  over  the  first  two  ion 
probes  (located  at  18.5  and  24.5  inches  respectively),  the  wave  speed  was  approximately 
1900  m/s  for  a  fill  fraction  of  0.5,  indicative  of  detonation.  The  wave  speed  calculated 
between  the  second  pair  of  ion  probes,  located  at  24.5  and  33  inches  respectively  was 
consistently  above  1800  m/s.  The  drop  at  the  third  ion  probe  was  not  nearly  as  sharp  due 
to  the  fact  that  it  was  located  just  behind  the  s-curve  in  the  detonation  tube.  As  the  shock 
wave  traveled  around  the  curve  of  the  tube,  it  began  to  round  as  did  the  combustion  front 
traveling  behind  it.  As  it  turned  around  the  s-curve  the  combustion  front  and  the  shock 
began  to  decouple,  meaning  that  the  detonation  was  weaker.  As  the  shock  rounded,  the 
perceived  wave  speed  decreased  slightly  and  the  voltage  drop  became  distorted.  The 
voltage  drop  across  the  ion  probe  was  distorted  because  a  decreased  number  of  ionized 
particles  were  passing  over  the  probe.  This  slower  wave  speed  is  indicative  of  a  weaker 
detonation. 
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Figure  IV-15:  Ion  probe  voltage  traces  for  ten  Hz  operation  at  desired  fill  fraction  of 
0.5  indicating  detonation  waves  speeds  at  ~70%  of  tube  length 


The  phenomenon  of  tube  heating  on  increasing  actual  fill  fraction  has  been 
documented  in  previous  experiments  where  using  a  desired  fill  fraction  of  0.34, 
detonation  wave  speeds  at  70%  tube  length  and  combustion  events  at  80%  tube  length 
were  documented  (Paxson,  et  al.  2011).  From  this  data  it  is  reasonable  to  assume  that 
even  at  a  desired  fill  fraction  of  0.5  the  actual  fill  fraction  was  well  over  0.7  and  possibly 
into  the  turbine. 

In  this  experiment,  at  a  fill  fraction  of  0.5  the  ion  probes  indicated  detonations 
occurring  at  a  length  corresponding  to  a  fill  fraction  of  0.69.  The  fact  that  detonations 
were  still  occurring  this  far  down  the  tube  implies  the  likelihood  that  combustion  and 
possibly  even  detonation  was  occurring  in  the  turbocharger  and  possibly  even  the  heat 
exchanger  itself  Combustion  and  or  detonations  occurring  inside  and/or  past  the  turbine 
would  have  the  effect  of  artificially  increasing  the  measured  thermal  efficiency  if 
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temperature  upstream  of  the  turbine  were  measured  direetly.  The  artifieially  increased 
efficiency  is  due  to  the  combustion  taking  place  after  the  flow  enters  the  turbine  or  after  it 
exits.  In  this  scenario  it  would  appear  that  the  temperature  at  the  exit  of  the  turbine  is 
higher  than  if  combustion  were  complete  before  the  flow  entered  the  turbine.  Burning 
into  or  past  the  turbine  may  even  cause  the  temperature  at  the  exit  of  the  turbine  to  be 
higher  than  the  temperature  at  the  turbine  inlet,  creating  efficiencies  greater  than  100%. 
To  the  contrary,  such  a  situation  would  actually  greatly  reduce  the  efficiency  of  the 
turbocharger,  as  the  turbine  would  not  have  the  opportunity  to  extract  energy  from  the 
flow  as  it  bums  after  passing  through.  To  prevent  this  issue  future  experiments  should 
include  an  ion  probe  at  the  entrance  and  possibly  even  at  the  exit  of  the  turbine  to 
determine  if  combustion  is  still  occurring. 


IV.5  Phase  Change 

It  was  learned  over  the  course  of  experimentation  that  during  PDC  operation,  the 
water  in  the  PDC  exhaust  products  was  condensing  inside  the  heat  exchanger.  The  fact 
that  the  water  was  condensing  meant  that  a  phase  change  was  taking  place  inside  of  the 
heat  exchanger  that  was  previously  unaccounted  for.  As  the  water  vapor  in  the  exhaust 
gas  condenses  on  the  surface  of  the  heat  exchanger  it  releases  energy  into  the  surface 
which  is  transferred  into  the  coolant  in  the  heat  exchanger.  This  energy  released  needed 
to  be  accounted  for  in  the  calculation  of  the  inlet  temperature  or  the  condensation  needed 
to  be  avoided. 
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Another  side  effect  of  the  vapor  condensing  on  the  heat  exchanger  walls  was  a 
reduction  in  heat  transfer  coefficient.  The  water  acts  as  another  barrier  between  the  hot 
gas  and  the  liquid  coolant  that  serves  to  remove  heat,  much  like  fouling  in  a  heat 
exchanger,  this  is  due  to  the  relatively  low  conductivity  of  water  (-.6  W/m  K)  in 
comparison  to  aluminum  (-180  W/m  K).  The  fact  that  the  water  releases  a  significant 
amount  of  energy  (2257  kJ/kg)  as  it  changes  phase  has  the  effect  of  overshadowing  the 
reduced  convection  coefficient  due  to  the  fouling  effect. 

In  an  attempt  to  avoid  condensation  on  the  inner  wall  of  the  heat  exchanger,  the 
length  of  the  heat  exchanger  was  halved.  This  shorter  heat  exchanger  gave  the  exhaust 
gas  less  surface  area  to  contact  and  transfer  heat  over,  resulting  in  increased  average 
exhaust  gas  temperatures.  This  hotter  flow  caused  increased  temperatures  at  the  heat 
exchanger  wall.  The  assumption  was  that  the  wall  temperature  would  be  above  the 
boiling  temperature  of  water  (373  K/212  F)  so  that  when  the  water  vapor  made  contact 
with  the  wall  it  would  not  condense  there. 

After  multiple  runs  with  the  ten  foot  long  heat  exchanger  over  varying  set  points, 
no  water  was  observed  exiting  the  heat  exchanger  during  equilibrium.  To  further  justify 
that  water  was  not  collecting  inside  the  heat  exchanger  a  simplified  analysis  was 
conducted.  Using  the  knowledge  that: 


q  =  mCpAT 


It  is  possible  to  calculate  the  heat  removed  by  the  water.  This  number  was  calculated  at 
ten  Hz  for  fill  fractions  of  0.5  and  0.8  without  the  turbocharger.  Due  to  the  relatively 
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small  thickness  of  the  aluminum  wall  and  the  high  eonduetivity  of  the  aluminum,  the 
temperature  gradient  across  the  wall  will  be  very  small,  so  it  ean  be  assumed  that  the 
temperature  on  one  side  of  the  wall  will  be  approximately  the  same  as  the  temperature  on 
the  other  side  of  the  wall.  To  determine  the  change  in  temperature  from  the  gas  flow  to 
the  hot  wall  the  following  formula  was  manipulated: 

q  =  hAAT 

where  the  heat  transfer  was  removed  by  the  water.  The  following  table  summarizes  the 
results: 


Table  IV-2:  Heat  transfer  comparison  for  fill  fractions  of  0.5  and  0.8 


fill  fraetion 

0.5 

0.8 

Qwater  (W^) 

22172 

32120 

h  (W/m^  K) 

350 

350 

AT  (K) 

87 

125 

average  hot  gas  temperature 
(K) 

603 

791 

average  hot  wall 
temperature  (K) 

516 

665 

Both  of  the  resulting  average  hot  wall  temperatures  are  well  above  the  boiling  point  of 
water  (373  K). 

The  temperatures  measured  at  the  exit  of  the  heat  exehanger  via  a  J-type 
thermocouple  supported  this  argument.  At  the  exit  of  the  heat  exchanger,  with  the 
turbine  in  the  flow,  exhaust  gas  temperatures  ranged  from  394  K  to  573  K  and  pressures 
measured  at  the  exit  of  the  heat  exchanger  were  slightly  below  atmospheric  ranging  from 
13.6  psia  to  14.2  psia.  The  temperatures  and  pressures  reeorded  would  not  support 
condensation. 
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The  lowest  temperature  noted  was  394  K  at  14.1  psia  when  the  PDC  was  running 


at  10  Hz  and  0.5  fill  fraetion.  At  this  pressure  the  water  in  the  exhaust  produets  would 
need  to  drop  to  371  K  or  23  K  below  the  temperature  observed  in  order  to  eondense.  Due 
to  the  temperature  differenee  from  the  gas  temperature  measured  in  the  centerline  of  the 
flow  to  the  wall  of  the  heat  exchanger,  there  may  have  been  water  condensing  inside  the 
heat  exchanger,  but  the  temperatures  at  the  exit  of  the  heat  exchanger  would  still  have 
been  high  enough  to  re-evaporate  any  water  that  might  have  condensed  on  the  walls.  Due 
to  the  control  volume  approach  of  this  technique,  the  water  condensing  and  evaporating 
inside  of  the  heat  exchanger  should  not  affect  the  calculated  temperature  at  the  inlet  of 
the  heat  exchanger. 


IV.6  Startup/Shutdown  Transient 

Another  phenomenon  that  was  exposed  during  experimentation  was  a  significant 
increase  in  heat  exchanger  exhaust  gas  temperature  after  detonations  were  terminated. 
One  example  of  this  occurrence  can  be  seen  in  Figure  IV- 16.  The  temperature  rose  for 
25  seconds  after  shutdown  before  it  began  to  cool.  Heat  exchanger  exhaust  gas 
temperature  peaked  at  23 1  F,  about  52%  higher  than  the  maximum  temperature  measured 
while  detonations  were  occurring.  For  all  experiments  with  the  original  20’  long  heat 
exchanger,  the  PDC  was  allowed  to  run  for  more  than  seven  minutes,  to  allow  for  a 
steady  measurement  to  be  recorded.  Upon  shutdown  of  the  PDC,  cooling  air  was  still 
flowing  through  the  detonation  tubes  at  the  same  temperatures  and  pressures  as  it  would 
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if  the  PDC  were  operating  without  a  spark.  This  temperature  spike  was  prevalent  in 
every  run  10  and  15  Hz  run  at  a  fill  fraction  of  0.5. 


24  Nov  10, 15  Hz,  No  Turbo 


Figure  IV-16:  Shutdown  transient  temperature  spike 


Though  the  exhaust  temperature  of  the  heat  exchanger  was  increasing  after 
detonations  ceased,  it  should  be  noted  that  the  water  temperature  did  decline  as  expected, 
and  the  overall  calculated  inlet  temperature  also  recedes.  The  calculated  inlet 
temperature  during  the  startup  and  shutdown  of  the  PDC  are  much  less  reliable  as  the 
experiment  has  not  reached  a  steady  state  condition,  therefore  the  energy  storage  term  is 
changing.  However,  it  is  not  the  value,  but  the  trend  that  is  important.  The  decline  in 
heat  exchanger  inlet  temperature  indicates  an  overall  reduction  in  the  energy  into  the 
system. 
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The  higher  temperature  readings  at  the  exhaust  of  the  heat  exehanger  after 
shutdown  are  unexpected  because  the  flow  entering  the  detonation  tube  is  slightly  above 
ambient  conditions  and  moving  at  subsonic  speeds.  The  slower,  colder  flow,  that  spends 
more  time  in  contact  with  the  heat  exchanger,  would  be  expected  to  produce  lower 
exhaust  temperatures  than  during  PDC  operation.  This  phenomenon  is  possibly  caused 
by  a  decreased  convective  heat  transfer  coefficient  after  shutdown.  Convective  heat 
transfer  coefficient  is  a  function  of  Nusselt  number,  which  is  lower  during  purge  and  after 
shutdown  (similar  to  purge)  than  during  the  blow  down  phase. 

There  are  a  number  of  factors  that  need  to  be  considered  with  the  convective  heat 
transfer  coefficient.  The  Nusselt  number,  seen  in  Eq.  5.1.9,  is  a  function  of  Reynolds 
number,  Prandtl  number  and  friction  factor.  The  Reynolds  number,  which  is  the  ratio  of 
inertial  forces  to  viscous  forces,  can  be  defined  as: 

Re  =  ^^^  Eq.  5.1.2 


In  the  transition  from  blow  down  to  purge,  L,  the  reference  length  (in  this  case  the  length 
of  the  wall  of  the  heat  exchanger)  does  not  change.  The  velocity  increases  by  an  order  of 
magnitude  and  the  density  decreases  slightly.  The  following  table  illustrates  a  simple 
comparison  of  the  two  conditions  to  show  why  the  heat  transfer  coefficient,  and  therefore 
the  energy  transferred  to  the  heat  exchanger  would  be  increased  during  blow  down: 
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Table  IV-3:  Simplified  comparison  of  blow  down  phase  and  purge  phase  heat 

transfer 


blow  down 

purge 

pressure  (Pa) 

506625 

101325 

temperature  (K) 

1200 

560 

density  (kg/m''3) 

1.47154932 

0.630664 

velocity  (m/s) 

30 

5 

conductivity  (W/m  K) 

0.078889806 

0.04389 

viscosity  (N  s/m''2) 

4.79474E-05 

2.9E-05 

Prandtl  number 

0.721889856 

0.689562 

Reynolds  Number 

1.40E+06 

1.66E+05 

friction  factor 

0.01098332 

0.016213 

Nusselt  number 

1530.321435 

262.9776 

convection  heat  transfer  coefficient  (W/m  K) 

1584.3407 

151.4709 

average  convection  heat  transfer  coefficient  (W/m  K) 

356.1665726 

heat  transferred  (W) 

421950.2515 

-30393.6 

total  heat  transfer  (W) 

34226.97906 

actual  heat  transferred  (from  water)  (W) 

-34000 

This  comparison  simplifies  the  analysis  by  assuming  a  uniform,  average 
temperature  of  the  purge  and  blow  down  gas  temperatures  and  the  wall  temperature.  The 
wall  temperature  is  assumed  to  be  835  K.  The  diameter,  circumference  and  surface  area 
used  in  this  calculation  were  all  based  off  of  the  ten  foot  heat  exchanger  used  in  later 
experimentation. 

In  this  simplified  comparison  the  pressure  during  blow  down  is  assumed  to  be 
approximately  five  atmospheres,  which  is  consistent  with  previous  results  (B.  H.  Schauer 
et  al.  2003),  while  the  purge  pressure  is  estimated  to  be  atmospheric.  The  temperature 
during  blow  down  is  estimated  to  be  1200  K  (a  conservative  estimate)  when  it  reaches  the 
heat  exchanger,  while  the  purge  temperature  used  is  the  same  as  what  is  calculated  in 
Appendix  C.  The  density  is  calculated  directly  assuming  the  working  fluid  is  an  ideal  gas 
(air)  with  gas  constant  286.9  (J/kg  K).  The  velocity  during  blow  down  is  estimated  to  be 
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six  times  faster  than  the  purge  flow.  The  conductivity,  viscosity  and  Prandtl  number  are 
all  based  off  of  the  associated  flow  temperature.  The  Reynolds  number,  friction  factor, 
Nusselt  number  and  heat  transfer  numbers  are  all  calculated  based  off  of  the  previously 
mentioned  values  and  the  equations  used  in  Appendix  A.  The  average  heat  transfer 
coefficient  and  average  heat  transferred  values  are  velocity  averaged,  meaning  that 
because  the  velocity  of  the  blow  down  is  six  times  higher  than  the  purge  velocity,  the 
blow  down  heat  transfer  coefficient  is  multiplied  by  one  seventh  and  the  added  to  the 
purge  heat  transfer  coefficient  which  is  multiplied  by  six  sevenths.  The  total  heat 
transferred  is  calculated  in  the  same  manner,  using  the  heat  transferred  during  each  phase 
and  multiplying  by  the  velocity  coefficient. 

The  resulting  heat  transfer  coefficient  and  heat  transfer  is  significantly  increased 
during  the  blow  down  phase.  Based  off  of  this  analysis,  the  heat  transferred  during  the 
purge  phase  is  actually  negative,  meaning  that  the  purge  air  is  cooling  the  inner  wall  of 
the  heat  exchanger.  It  should  be  noted  that  the  average  heat  transferred  (34.2  kW)  is  in 
the  range  of  what  was  calculated  to  be  removed  by  the  water  during  testing  which  was 
between  22  kW  and  35  kW. 

It  is  also  hypothesized  that  the  pressure  oscillation  had  a  great  deal  to  do  with  the 
temperature  spike  at  shutdown.  The  temperature  spike  was  more  prevalent  with  the 
twenty  foot  long  heat  exchanger.  The  pressure  oscillation  was  constantly  pulling  cool  air 
back  into  the  heat  exchanger  for  a  larger  portion  of  every  cycle  during  operation,  creating 
a  cooling  effect  over  the  thermocouple.  When  the  PDC  stopped  detonating,  the  pressure 
oscillation  almost  completely  disappeared  and  as  a  result  the  thermocouple  is  in  contact 
with  a  more  consistent  flow  temperature  associated  with  the  purge  flow  that  is  heated  by 
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the  detonation  tube.  This  flow  temperature  is  aetually  hotter  than  what  the  thermocouple 
at  the  heat  exchanger  exit  was  in  contact  with  on  average  during  operation. 

The  temperature  spike  was  significantly  reduced  and  observed  when  operating 
with  the  9.5  foot  extension  at  the  exit  of  the  heat  exchanger  at  lower  frequencies  and  fill 
fractions.  In  the  same  configuration  at  higher  frequencies  and  fill  fractions  the 
temperature  spike  was  not  observed  was  not  noticeable.  The  temperature  spike  was  not 
present  because  the  average  flow  temperature  registered  by  the  thermocouple  was  greater 
than  the  purge  flow  temperature,  even  after  PDC  tube  heating. 


IV.7  Analytical  Errors  and  Corrections 

After  stepping  back  from  analysis,  it  appears  that  there  were  four  main  causes  of 
error  in  the  analysis  of  data  for  this  experiment:  pressure  oscillations  and  flow  speeds, 
water  condensing  in  the  heat  exchanger,  water  mass  flow  rate  inaccuracy  and  water  outlet 
temperature  precision.  These  issues  were  evaluated  at  10  Hz  and  a  fill  and  purge  fraction 
of  0.5  and  0.9  because  this  is  where  the  largest  disagreement  was  found  between  the 
calculated  and  measured  values.  The  evaluation  looks  at  a  maximum  value  that  would  be 
expected  if  the  error  were  corrected  and  possible  solutions  to  correct  for  the  error,  either 
analytically  or  experimentally.  The  results  of  this  evaluation  are  located  in  the  table 
below: 
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Table  IV-4:  Leading  causes  of  error,  adjusted  values  and  possible  solutions 


issue 

temperature  (K) 

change 

possible  solution 

thermocouple  contact 
ratio 

932 

(+39%  from  670  K) 
(+4.5%  from  892  K) 

1:1  vs.  1:6 
blow  down:  purge 

increase  frequency/mass 
flow  rate, 

heat  exchanger 
condensation 

553 

(-18%  from  670  K) 

0%  to  100%  relative 
humidity 

increase  exhaust  gas 
temperature 

water  flow  rate 

688  K  to  664  K 
(+2.7%  to  -0.9%  from 
670  K) 

+.25  gal/min  to  -.25 
gal/min 

increase  water  flow  rate, 
increase  accuracy  of 
meter 

water  temperature 

738  K  to  614  K 
(+10%  to  -8.4%  from 
670  K) 

+2  K  to  -2  K 

decrease  water  flow  rate 

Original  calculated  heat  exchanger  inlet  temperature  was  670  K  and  the  measured  value  was  892  K 

The  first  and  largest  cause  of  error  evaluated  was  the  lack  of  accounting  for  the 
variance  in  the  speed  of  the  flow  and  the  ratio  of  time  each  portion  of  the  flow  spends  in 
contact  with  the  thermocouple.  Due  to  the  determination  that  the  flow  is  still  segregated 
as  it  passes  over  the  thermocouple,  the  blow  down  and  purge  flows  spend  a  varying 
length  of  time  in  contact  with  the  thermocouple.  Through  observation  of  the  Schlieren 
video  it  was  determined  that  the  blow  down  gas  spends  approximately  1/6*  of  the  cycle 
in  the  frame,  while  the  rest  is  significantly  cooler  gases.  It  is  reasonable  to  assume  that 
the  flow  at  the  exit  of  the  heat  exchanger  exhaust  pipe  is  very  similar  to  what  is  seen  at 
the  thermocouple  at  the  exit  of  the  heat  exchanger,  therefore  this  was  the  assumption.  A 
very  simplified  analysis  was  used  to  determine  a  corrected  value.  The  evaluation  started 
by  estimating  the  temperature  of  the  blow  down  gas  to  be  2500  K  based  off  of  recent  high 
speed  measurements  (K.  P.  Rouser,  P.  I.  King  and  F.  I.  Schauer,  et  al.  2011).  It  was  also 
estimated  that  the  purge  gas  entered  the  heat  exchanger  at  the  temperature  approximated 
by  the  original  calculations  of  382  K.  At  this  point  it  was  estimated  that  all  of  the  energy 
transferred  into  the  water  was  done  so  by  the  blow  down  gas.  The  energy  transfer  would 
bring  the  temperature  of  the  blow  down  gas  to  1945  K.  The  purge  temperature  would 
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likely  increase  slightly,  removing  energy  from  the  heat  exchanger  inner  wall,  so  it  is 
estimated  to  be  400  K  when  it  reaches  the  thermocouple.  Applying  the  ratio  previously 
mentioned  it  is  possible  to  calculate  an  average  heat  exchanger  exit  temperature  yields 
657  K  versus  what  was  measured  at  393  K: 


T  =(T  if— 1  +  fr  'if— 

avg  V  blowdown  J  /^  \  \  purge  J  ^ 

Hx  \oJ  yoj 


=  (1945  K)  -  +(400  .S:)  -  =651  K 

V  6  y  V  6  ; 


This  temperature  is  now  applied  to  the  original  heat  exchanger  analysis  via  the 
method  described  in  section  IIL4.  The  resulting  average  heat  exchanger  inlet  temperature 
after  applying  this  correction  is  932  K.  This  value  is  approximately  4.5  %  higher  than  the 
measured  temperature,  and  is  in  much  better  agreement  with  the  measured  value. 

This  issue  could  be  corrected  analytically  as  discussed  in  the  previous  paragraph  or 
experimentally  by  increasing  the  frequency  of  operation,  and/or  mass  flow  rate  of  the 
fuel-air  mixture.  Increasing  frequency  would  increase  the  speed  that  the  purge  gas  is 
pushed  through  the  tube  by  the  blow  down,  and  the  ratio  would  be  closer  to  1:1. 
Increasing  the  mass  flow  rate  of  the  main  fiiel-air  mixture  would  increase  the  volume  of 
the  blow  down  phase  and  as  a  result  the  thermocouple  would  be  in  contact  with  the  hot 
gas  for  a  longer  period  of  time  and  record  an  increased  temperature.  The  main  fdl 
fraction  was  varied  during  the  experiment  and  this  increasing  temperature  was  verified. 

The  next  issue  that  was  analyzed  was  the  possibility  of  water  condensing  in  the  heat 
exchanger.  In  order  to  approximate  the  difference  in  average  heat  exchanger  inlet 
temperature  if  condensation  were  occurring  the  worst  case  scenario  was  considered  if 
100%  relative  humidity  were  measured  at  the  exit  of  the  heat  exchanger.  The  method  for 
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calculating  this  temperature  is  described  thoroughly  in  Appendix  E.  The  resulting 
temperature  if  all  of  the  water  in  the  exhaust  products  condensed  into  liquid  inside  of  the 
heat  exchanger  was  533  K  versus  the  670  K  measured.  To  correct  the  issue  of  water 
condensing  in  the  heat  exchanger,  the  exhaust  gas  temperature  would  need  to  increased 
so  that  the  water  in  the  products  remains  well  above  373  K  (212  F),  the  temperature 
where  the  water  would  condense.  The  experiment  was  actually  modified  to  correct  this 
issue  and  the  resulting  (uncorrected)  heat  exchanger  exit  temperatures  were  above  370  K. 

To  approximate  the  error  caused  by  the  inaccuracy  of  the  mass  flow  rate 
measurement,  the  original  analysis  from  Section  III.4  was  conducted,  but  water  mass 
flow  rate  was  increased  and  decreased  by  .25  gallons  per  minute  (the  precision  error  of 
the  meter).  The  increase  in  water  flow  rate  yielded  an  average  heat  exchanger  inlet 
temperature  of  688  K  and  decreasing  it  yielded  a  temperature  of  664  K.  To  reduce  this 
error  the  mass  flow  rate  of  the  water  could  be  increased.  Increasing  the  water  mass  flow 
rate  would  serve  to  reduce  the  error  in  proportion  to  the  measured  rate.  Increasing  the 
mass  flow  rate  would  also  have  the  negative  side  effect  of  increasing  the  percentage  of 
precision  error  of  water  outlet  temperature  and  as  a  result,  it  is  recommended  that  a  more 
accurate  water  flow  meter  be  acquired. 

Finally  to  quantify  the  precision  error  of  the  water  outlet  temperature  a  similar 
analysis  as  discussed  in  the  previous  paragraph  was  conducted.  Increasing  the  water 
temperature  by  two  Kelvin  increased  the  calculated  average  heat  exchanger  inlet 
temperature  to  738  K  and  the  opposite  reduced  the  temperature  calculated  to  614  K.  To 
reduce  the  percentage  of  the  precision  error  on  water  outlet  temperature  the  water  flow 
rate  could  be  reduced.  Reducing  the  water  mass  flow  rate  would  raise  the  change  in 
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water  temperature  from  the  inlet  to  the  exit  of  the  heat  exchanger  and  decrease,  by 
percentage,  the  precision  error. 
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V.  Conclusions 


One  of  the  original  assumptions  made  prior  to  this  experiment  was  that  the 
temperatures  and  pressures  calculated  at  the  inlet  of  the  heat  exchanger  when  the 
turbocharger  was  not  present  would  be  representative  of  the  turbine  inlet  conditions  when 
the  turbocharger  was  present.  This  analysis  also  hypothesized  that  the  calculated 
temperatures  and  pressures  at  the  inlet  of  the  heat  exchanger  would  be  representative  of 
the  turbine  exit  temperature  when  the  turbine  was  in  the  flow.  It  became  very  clear  upon 
examining  the  data  that  this  assumption  was  not  accurate  for  the  conditions  analyzed.  It 
was  found  that  the  pressure  oscillations  at  the  exit  of  the  heat  exchanger  were  more 
pronounced  without  the  turbocharger  to  act  as  a  damper.  The  increased  pressure 
oscillations  caused  for  even  lower  temperature  measurements  when  the  turbocharger  was 
not  present.  The  inaccuracy  was  so  extreme  in  fact  that  the  results  produced  turbine  exit 
temperatures  greater  than  the  turbine  inlet  temperatures. 

The  addition  of  the  turbocharger  also  adds  back  pressure  to  the  PDC  tube  and  as  a 
result,  increased  pressures  and  temperatures  are  observed.  The  trend  of  increasing 
combustion  temperatures  with  increasing  pressure  is  consistent  with  analytical  results 
(Schultz  2000).  Manifold  pressures  increased  by  approximately  10%  with  the  addition  of 
the  turbocharger.  Combustion  occurring  in  and  through  the  turbine  could  also  have 
contributed  to  the  increased  temperatures  measured  at  the  exit  of  the  turbine  in 
comparison  to  without  the  turbine. 

Significant  pressure  oscillations  were  observed  during  experimentation  at  the  exit 
of  the  heat  exchanger  which  was  likely  the  cause  for  the  lack  of  conformity  between  the 
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calculated  temperatures  and  those  that  were  measured  directly.  It  is  likely  that  at  higher 
frequeneies,  approximately  65  Hz  and  greater,  the  pressure  oscillations  would  be 
attenuated  and  the  caleulated  inlet  temperature  of  the  heat  exehanger  would  be  in  eloser 
agreement  with  the  mass  averaged  temperatures  at  the  same  location.  Operation  at  higher 
frequeneies  is  also  supported  by  the  fact  that  at  10  and  15  Hz  operation,  the  turbine 
remains  unchoked  for  a  majority  of  the  cycle.  At  higher  frequeneies,  the  opposite  would 
be  true  and  likely  inerease  effieieney. 

Originally  the  test  matrix  for  this  experiment  included  operation  at  higher  and 
more  varied  frequencies  (10,  15  and  20  Hz).  Unfortunately  the  PDC  began  baekfiring 
during  extended  runs.  Backfiring  was  due  to  auto-ignition  where  the  spiral  used  for  DDT 
became  so  hot  that  the  fuel-air  mixture  began  igniting  before  the  valves  elosed  behind  it. 
This  diseovery  foreed  extended  runs  to  be  conducted  at  lower  frequencies.  It  was  noted 
that  the  tube  temperatures  remained  higher  when  the  turboeharger  was  applied  at  the  exit 
of  the  PDC  tube.  Use  of  alternate  fuels,  like  ethylene,  may  allow  the  PDC  to  be  run  for 
extended  periods  at  higher  frequeneies. 

Prior  to  this  experiment  it  was  hypothesized  that  the  addition  of  the  heat 
exchanger  and  turbocharger  would  serve  to  help  mix  and  steady  the  flow.  After 
analyzing  the  Schlieren  high  speed  videos  it  was  apparent  that  the  flow  was  still 
extremely  unsteady  at  the  exit  of  the  heat  exehanger.  The  combination  of  increasing 
frequency  and  adding  the  turbocharger  resulted  in  reduced  pressure  oscillations  at  the  exit 
of  the  heat  exehanger,  but  the  flow  still  pulsed  into  and  out  of  the  exhaust  pipe  several 
times  over  the  course  of  one  cycle  at  20  Hz. 
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Though  uncertainty  values  for  turbine  work,  compressor  work  and  efficiency 
were  higher  than  desired,  the  resulting  values  are  still  reasonable.  For  a  turbine  powered 
by  unsteady  flow  that  is  unchoked,  efficiencies  between  10  to  15%  are  not  outlandish. 
Ideally  a  time  accurate  efficiency  could  be  attained  and  integrated  over  an  entire  cycle  to 
provide  a  comparison. 

The  process  of  using  compressor  mass  flow  rates,  temperatures  and  pressures  to 
estimate  compressor  work  was  found  to  be  consistent.  The  calculated  turbine  exit 
temperatures  agree  more  with  those  that  were  measured  directly  as  frequency  increased. 
The  temperatures  measured  using  thermocouples  at  the  exit  of  the  PDC  and  turbine  were 
closer  to  the  mass  weighted  average  temperatures  than  those  calculated  however,  as  PDC 
frequency  increases  a  thermocouple  is  not  likely  to  survive,  let  alone  measure 
temperature  accurately.  The  energy  balance  approach  would  be  best  suited  for  a  scenario 
where  PDC  frequency  is  higher,  in  order  to  reduce  the  likelihood  that  pressure 
oscillations  allow  ambient  air  to  reach  the  thermocouple  at  the  heat  exchanger  exit. 

A  mass  averaged  temperature  is  necessary  for  calculating  work  and  efficiency. 
The  temperature  recorded  by  a  thermocouple  is  actually  the  result  of  heat  transfer  to  and 
from  the  thermocouple  during  unsteady  flow.  This  thermocouple  measured  temperature 
will  fall  short  of  the  mass  averaged  temperature  for  an  unsteady  flow. 
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VI.  Recommendations 


Future  experimentation  should  include  operation  at  higher  frequencies  and  mass 
flow  rates  for  reasons  mentioned  previously.  The  majority  of  the  error  in  calculated 
turbine  inlet  and  exit  temperatures  was  the  result  of  pressure  fluctuations  at  the  exit  of  the 
heat  exchanger.  More  emphasis  should  be  placed  on  steadying  the  flow  of  the  PDC  to 
acquire  a  more  accurate  measurement.  This  could  be  accomplished  by  using  a  larger 
vessel  (on  the  order  of  ten  times  the  volume  of  the  heat  exchanger  used  in  this 
experiment)  at  the  exit  of  the  turbine  where  the  shock  wave  would  have  room  to  expand 
and  dissipate.  Without  steadying  the  flow  there  is  little  hope  for  measuring  a  mass 
averaged  temperature.  Schlieren  or  other  flow  visualization  techniques  should  be  used  to 
verify  steady  flow  at  the  exit  of  the  device.  Arrangements  should  be  made  to  account  for 
any  burning  in  or  through  the  turbine,  especially  during  extended  runs.  Using  ion  probes 
immediately  in  front  of  and  behind  the  turbine  would  help  to  identify  any  combustion  that 
may  be  occurring  in  the  area.  Any  future  approaches  to  measuring  enthalpy  at  the  inlet  or 
exit  of  the  turbocharger  would  need  to  steady  the  flow  before  a  temperature  measurement 
is  taken. 
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VII.  Appendices 


Appendix  A  -  Heat  Exchanger  Design 

The  following  formulas  and  methods  were  used  in  combination  with  Microsoft  Excel®  in 
order  to  produce  estimates  of  heat  exchanger  parameters  (length,  diameter,  etc.). 

First,  an  initial  guess  at  inlet  temperature  and  desired  outlet  temperature  for  PDC 
the  hot  gases  were  necessary.  Initial  conservative  estimates  for  heat  exchanger  inlet 
temperature  of  1450  K  (2150  F)  and  outlet  temperature  of  900  K  (1 160  F)  were  used. 

This  inlet  temperature  was  greater  than  what  the  expected  average  inlet  temperature 
would  be  and  the  outlet  temperature  was  selected  to  ensure  the  exhaust  gases  would  be 
cold  enough  to  be  directly  measured  with  a  J-type  thermocouple  (<1382  F).  A  total  mass 
flow  of  5.3  Ib/min  was  selected  and  the  constant  pressure  specific  heat  for  the  inlet  and 
outlet  flow  is  attained  from  a  source  (Incropera,  et  al.  2007).  From  these  numbers,  q  is 
calculated: 

q  =  Thg.out)  Eq.  5.1.3 

All  of  the  heat  is  assumed  to  go  straight  into  the  water: 

qhg=qwater=q  Eq.  5.1.4 

The  outlet  temperature  of  the  water  is  calculated  in  the  reverse  manner,  based  off  of  the 
water  mass  flow  rate  and  the  water  inlet  and  outlet  temperatures  and  specific  heats.  The 
open  loop  cooling  system  has  a  constant  liquid  inlet  temperature.  The  mass  flow  rate  of 
the  liquid  is  manipulated  to  maintain  an  outlet  temperature  that  is  below  boiling: 
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+  C  T 

Pliq,in 


^liq,out 


C 


Eq.  5.1.5 


Pliq  ,01 


Next  the  Reynolds  numbers  of  the  hot  gas  flow  and  the  liquid  flow  must  be  determined: 

4m, 


''hg 


Eq.  5.1.6 


and 


Re,,  = 


^  (Do +D,- )/«/,■, 


Eq.  5.1.7 


Because  the  coolant  flows  through  the  annulus,  the  same  calculation  is  not  used  for 
Reynolds  number.  The  friction  factor  was  calculated  using: 


/  =  (0.790  In  (Re) -1.64)“"  Eq.  5.1.8 

The  hot  gas  Reynolds  number  is  above  10,000  so  the  correlation  for  Nusselt  number 
(Nu),  the  ratio  of  convection  to  pure  conduction  heat  transfer,  used  is: 


(//8)(Re-1000)Pr 
l  +  12.7(//8)"^(Pr^^^-l) 


Eq.  5.1.9 


Prandtl  number,  a  ratio  of  the  momentum  diffusivity  (or  v  the  kinematic  viscosity  which 
is  the  dynamic  viscosity  divided  by  density)  to  thermal  diffusivity  (a),  is  assumed  to  be 
that  of  air.  The  friction  factor  and  Nusselt  number  of  the  liquid  is  calculated  in  the  same 
fashion.  The  convective  heat  transfer  coefficient  for  the  hot  gas  side  is  then: 

Eq.  5.1.10 


and  the  convective  heat  transfer  coefficient  for  the  liquid  side  is: 
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Next  the  overall  heat  transfer  coefficient  is  calculated: 


U  = 


J-+-L 

hfe  hi,, 


The  next  step  is  to  calculate  the  log  mean  temperature  difference: 


AT.„= 


(t  _T  "l-fx  _T  ] 

_  \  hg,in  liq,in  J  \  hg,out  liq,out  / 


[("^hg,in  ■Tliq.in  )  ^  (Thg^out  "Tuq.out  )] 


Eq.  5.1.11 


Eq.  5.1.12 


Eq.  5.1.13 


Now  all  of  the  necessary  values  have  been  attained  to  allow  for  the  calculation  of  the 
overall  length  of  the  system: 


L  = 


U^D,AT,„ 


Eq.  5.1.14 


A  snapshot  of  the  spreadsheet  used  for  this  analysis  is  seen  in  Figure  VII- 1 : 


Given  SI  Units 

Concentric  tube  counterflow 

US  units 

Input  Variables/Fixed  Values 
Table  Lookup 

Output  Variables 

1 

T,hg,in  (K)  14S0 

2150.6  deg  F 

T,hg,out  (K)  900 

1160.6'  deg  F 

mdot*Cciout*Tout-mdc: 

1 

mdothg  (Kg/s)  0.04 

5.291066667 

Ib/min 

q,hg(W) 

Re,D 

Nu,D 

h,hg 

Cp,in,hg(J/Kg*K)(air)  1225 

Pr.avg 

0.685 

30694 

11015.2549 

31.86145992 

32.4879493 

Cp,out.hg  (J/Kg*K)  (air)  1121 

form  8.62  (Incropera) 

friction  factor 

Cf 

T,iiq,in  (K)  293 

68 

degF 

0.030643519 

0.00766088 

TJiq^out  (K)  313.75 

105.3500615' 

degF 

Cp,in,liq  (J/Kg*K)  4181 

u,,!iq=q,h? 

Cp,oirt,liq  (J/Kg*K)  4184 

q,liq(W) 

liquid  Re,D 

Nu,D 

h,liq  (W/m*2K) 

mdotjiq  (Kg/s)  0.35 

5.548216645 

gal/min 

30694 

2.73E+03 

1.71E+01 

1.51E402 

: i*Dhvd*2  4 

1 

D,tube,in  (m)  0.089 

3.50393712 

inch 

deltaTJm  (K) 

Ac  (m*2) 

Ac(in''2) 

Length  (in) 

240 

D,tube,out  (m)  0.1615 

6.35826792 

inch 

844.152919 

0.014263812 

22.10895485 

Liauid  Pronerties 

D,hyd  (m)(iiquid)  0.0725 

2.8543308 

inch 

visc,avg  (Ns/mA2)  (iiqui  0.000651 

Vise  in 

vise  out 

Prin 

Prout 

Kin 

Kout 

Pr,avg  (liq)  4.38 

9.59E-04 

3.43E-04 

6.62 

2.14 

0.606 

0.671 

k,avg  (W/mK)  (liq)  0.6385 

Hot  Gas  Properties 

visc,avg  (Ns/mA2)  5.2E-05| 

Vise  in 

vise  out 

Prin 

Prout 

Kin 

Kout 

PrJH 

5.90E-05 

4.49E-05 

0.687 

0.728 

0.11 

0.0715 

k,avg  (W/mK)  0.09075| 

Process 

U  (Overall  Heat  Transfer  Coefficient)  (W/m''2K) 

2.67E40i| 

friction  factor 

4.70E-02 

Cf 

1.18E-02 

Mdot*Cp*{Tin-Tout)=qi  qhg=-qljq 

qliq=U*A*Tim 

A=pi*Din‘L 

L,  pipe  length  (m) 

in  feet 

1 

Dhyd=Dout-Din 

5.97E400 

1.96E401 

1  1 

Figure  VII-1:  Sample  heat  exchanger  design  spreadsheet 
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It  is  notable  from  Figure  VII- 1  that  the  convective  heat  transfer  coefficient  for  the  liquid 
side  is  much  larger  than  for  the  hot  gas  side.  This  translates  physically  into  the  wall  of 
the  heat  exchanger  being  primarily  the  same  temperature  as  the  liquid  flowing  past  the 
wall. 


Table  VII-1:  Original  heat  exchanger  properties 


Annulus  Volume 

.0895  m^ 

3.16  ft^ 

Inner  diameter  volume 

.0328  m^ 

1.159  fl^ 

Volumetric  gas  flow  rate 

845.7  mVs 

147.2  ftVs 

Volumetric  coolant  flow 
rate 

3.47  mVs 

.735  flVs 

Average  coolant  velocity 

5.411  E'^m/s 

1.775  E'^ft/s 

Average  air  velocity  (10 
Hz) 

13  m/s 

9.845  fl/s 

Coolant  cycle  time 

1126  s 

18.77  min 

Hot  gas  cycle  time 

0.5  s 

.03  min 
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Appendix  B  -  Free  convection  and  radiation  calculations  for  PDC  tube 


and  heat  exchanger 


The  following  formulas  and  methods  (Incropera,  et  al.  2007)  were  used  in 
combination  with  Microsoft  Excel®  in  order  to  produce  estimates  of  heat  exchanger 
losses  to  include  free  convection  and  radiation. 

The  heat  transfer  at  the  boundary  of  the  system  (out  of  the  heat  exchanger)  is: 

=ii^D(T.-Tj  +  CTD<T(T.‘-T„‘)  Eq.  5.1.15 


During  operation  the  highest  temperature  the  surface  of  the  heat  exchanger  reached  was 
105  F  (313  K)  so  a  conservative  estimate  of  120  F  (322  K)  is  used  for  Ts,  surface 
temperature,  in  this  analysis.  To  calculate  the  heat  transfer  coefficient  for  free  convection 

{h),  Rayleigh  number  (Ra),  the  ratio  of  the  buoyancy  forces  to  the  viscous  forces  in  the 
fluid,  must  first  be  calculated: 


Ra 


gA(T.-T.)D^ 

va 


Eq.  5.1.16 


(Note:  in  the  case  of  free  convection  air  is  assumed  to  be  an  ideal  gas  and  an  average 
value  of  g,  p,  V  and  a  is  used.)  For  an  ideal  gas: 

1 


Eq.  5.1.17 


avg 


The  empirical  correlation  for  Nusselt  number  of  a  long  horizontal  cylinder  (which  the 
heat  exchanger  will  be  approximated  as)  is: 
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2 


0.60 +  - 


0.387i?a 


1/6 


1  +  (0.559 /Pr)' 


9/16 


8/27 


Eq.  5.1.18 


Finally  free  convection  heat  transfer  coefficient  is  calculated  in  the  same  manner  as  in 
Eq.  5.1.10. 

The  numbers  for  the  radiation  calculation  qrad  were  directly  inserted  into 
Eq.  5.1.15.  It  should  be  noted  that  qtotai  is  the  sum  of  qrad  and  qnat  conv  multiplied  by  the 


length  of  the  heat  exchanger. 


Heat  Transfer  from  Heat  X 

F 

K  ^ 

V‘ 

gas) 

g(m/s''2) 

sigma  (W/m''2  KM) 

9.8 

5.67E-08 

Ttube,avg  (conservative)  123 

323.705556 

Tamb  65 

291.483333  " 

ua 

Tdelta  58 

32.2222222 

Tavg 

307.594444  . 

p^n^  0.387i?a,” 

K(.559'Prf” 

rho,air,avg  (kg/m''3) 

1.1614  u 

1 

Beta  (K^'-l) 

0.03103448 

D(ln,m)  6.625* 

Utube(in,m)  240 

0.168275  1  ^ 

,+9',^  = 

kDil- 

■L 

■C) 

nu,avg(m^2/s) 

alpha,avg(m''2/s) 

Pr,avg 

k,avg  (W/mK) 
emiss,AI 

1.59E-05 

2.25E-05 

0.707 

2.63E-02 

0.2 

Ra,0 

1.31E+08 

Nu,D 

61.3771738 

h(W/m''2K) 

9.592748 

q',nat.conv.  (W/m) 

163.405999 

q',rad.  (W/m) 

22.5487927 

q,total  (Hp,  W)  1.518998 

1133.58041 

mdot  (air) 

Ib/min 

kg/s 

7.84 

0.059192 

Figure  VII-2:  Sample  free  convection  and  radiation  calculation  spreadsheet 
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Appendix  C  —  Purge  gas  heat  exchanger  inlet  temperature 


As  the  purge  gas  travels  through  the  PDC  detonation  tube,  heat  is  transferred  from 
the  tube  to  the  gas  in  a  cooling  process.  This  change  in  temperature  is  significant  and 
cannot  be  ignored.  The  following  process  provides  an  estimate  of  the  purge  gas 
temperature  as  it  enters  the  heat  exchanger. 

First  the  manifold  temperature  (Tman)  is  measured  directly.  This  is  the 
temperature  of  the  air  as  it  enters  the  PDC  tube.  The  PDC  tube  temperature  is  also 
measured  at  the  front  and  back  ends  of  the  tube  to  determine  an  average  tube  surface 
temperature  (Tsurf).  The  mass  flow  rate  (m)  of  the  purge  flow  is  directly  measured  as 
well.  The  Reynolds  number  and  friction  factors  are  determined  using  Eq.  5.1.6  and 
Eq.  5.1.8.  The  Nusselt  number  is  found  by  again  using  Eq.  5.1.9  for  fully  developed 
turbulent  flow  in  a  smooth  circular  tube.  Next  the  convective  heat  transfer  coefficient  is 
calculated  with  Eq.  5.1.10.  Finally  the  purge  gas  temperature  at  the  heat  exchanger  inlet 
is  calculated,  assuming  a  constant  surface  temperature  based  off  of  the  average  tube 
surface  temperature  (Tsurf)  and  : 


in 


Eq.  5.1.19 


All  properties  of  air  are  based  off  of  curve  fits  from  tabular  data  (Incropera,  et  al.  2007), 
to  include  specific  heat,  thermal  conductivity,  Prandtl  number  and  viscosity.  An  average 
temperature  of  the  purge  gas  is  used  to  find  the  air  properties.  For  the  reason  that  the 
values  of  the  properties  vary  with  temperature,  this  is  an  iterative  process. 
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T,man/T,in  (K)  | 

1  2931 

T,surf  (K) 

1200^ 

m,dot,purge 

0.015787 

vlsc,purge  (N  s/m''2) 

2.39885E-05 

Tube  Length  (m) 

1.2192 

Tube  Diam  (m) 

0.0547878 

Re 

f 

15294.0818" 

0.028040509 

Pr 

0.690646642 

n  (.4  heating,. 3  cooling) 
Nu 

0.4 

44.16201946  41.41097 

k(W/mK) 
h  (W/(m''2K) 

Cpiair) 

0.036562142 

27.63523417 

1017.765672 

lr,DurRe,out(K| 

567.790Hii 

T,avg,purge  (K) 

430.42  430.3983 

Figure  VII-3:  Sample  purge  gas  inlet  temperature  calculation  spreadsheet 


86 


Appendix  D  —  Original  Experiment  and  Modifications 


Due  to  the  many  lessons  learned  in  the  process  of  conducting  this  experiment,  it 
was  necessary  to  continually  adapt  and  modify  the  experiment  to  attain  accurate  turbine 
inlet  temperatures.  This  appendix  includes  a  summary  of  those  changes. 

Originally  two  aluminum  pipes  were  cut  from  single  20  foot  sections  into  two  10  foot 
sections  for  ease  of  transportation  and  to  prevent  the  experiment  from  protruding  over  the 
side  of  the  test  stand.  A  %  inch  diameter  bent  tube  allows  water  to  flow  between  ten  foot 
sections  of  the  heat  exchanger.  A  180°  mild  steel  elbow  turned  the  exhaust  products  back 
into  the  second  section  of  the  heat  exchanger.  At  this  point  the  exhaust  products  are 
aimed  back  at  the  PDC  rig,  so  a  90°  steel  elbow  turns  the  flow  again  where  the 
temperature  of  the  gas  was  measured  eight  inches  downstream  of  the  heat  exchanger  exit. 
A  simplified  schematic  of  the  heat  exchanger  from  a  top  view  is  pictured  in  Figure  VII-4: 


10' 


H20in.  -70F, 


-flZ 

TurboiPDE 

exIiaiJGt,  -^200DF 

tH20 
out.  3/4=’ 


_  -'S.Bgal/min 
I  j/4  ’  inlet 


6’  Diain,  Soiled  1  [) 
606 1T6 


exit 


Figure  VII-4:  Top  view  of  heat  exchanger 
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Figure  Vll-5  shows  the  flow  direction  for  both  hot  gas  and  coolant  through  the  heat 


exchanger  and  Figure  Vll-6  identifies  the  locations  of  thermocouples. 


Figure  VII-5:  Counter  flow  heat  exchanger  attached  directly  to  PDC 
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Figure  VII-6:  Heat  exchanger  instrumentation 


The  eleven-blade,  radial  turbine  used  in  initial  experimentation  was  a  Garrett  T3 
automotive  turbocharger,  pictured  in  Figure  Vll-7.  The  turbine  wastegate  was  capped  so 
that  all  combustor  exhaust  passes  through  the  turbine.  The  T3  was  also  equipped  with  a 
radial  compressor  having  six  primary  impeller  blades  and  six  splitter  blades.  The  water- 
cooled,  center  housing  contains  the  shaft  and  dual  journal  bearing  assemblies.  The  T3  has 
a  0.58  A/R,  45  trim  compressor,  a  T3  4-Bolt  inlet  and  T3  5-bolt  discharge  exhaust. 
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Figure  VII-7:  Garrett  T3/T4E  automotive  turbocharger 


Figure  VII- 8  illustrates  the  T3  eompressor  operating  map  where  target  eompressor 
operating  eonditions  run  down  the  eenter  of  the  effieieney  islands. 
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Figure  VII-8:  Garrett  T3  compressor  operating  map 
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Originally  the  experiment  was  run  in  two  different  configurations:  configuration 
1,  without  a  turbocharger  seen  in  Figure  VII-9  and  configuration  2,  with  a  turbocharger 
seen  in  Figure  VII- 10.  Both  configurations  used  a  20  foot  long  heat  exchanger. 


Figure  VII-9:  Configuration  I-PDC  and  heat  exchanger 
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Figure  VII-10:  Configuration  2-PDC,  turbocharger  and  heat  exchanger 


Preliminary  results  showed  that  due  to  low  heat  exchanger  exhaust  temperatures  that 
the  water  in  the  exhaust  products  was  able  to  condense  inside  the  heat  exchanger.  To 
reduce  the  amount  of  energy  removed  from  the  exhaust  and  prevent  water  from 
condensing  the  heat  exchanger  length  was  reduced  to  ten  feet. 

The  majority  of  the  results  derived  in  this  paper  are  the  result  of  configuration  two, 
with  a  ten  foot  heat  exchanger.  It  was  determined  that  the  data  produced  without  the 
turbocharger  did  not  mimic  accurately  the  conditions  at  the  inlet  of  the  turbocharger. 
The  increased  pressure  oscillations  noted  upon  removing  the  turbocharger  we  the  most 
likely  cause  of  the  reduced  accuracy  of  this  approach. 


92 


Appendix  E—  Data  analysis  with  water  condensing 


During  the  experiment  using  a  20  foot  long  heat  exchanger  it  was  noted  that  water 
condenses  inside  the  heat  exchanger.  Condensing  water  signifies  a  phase  change  is 
taking  place.  The  phase  change  suggested  additional  energy  was  being  released  into  the 
wall  of  the  heat  exchanger  which  was  not  accounted  for  in  previous  calculations.  Starting 
from  Eq.  2.3.3,  a  latent  heat  term  must  be  added  to  account  for  the  phase  change: 


E  E  —  E  E  “h  E 

coolant  gas  coolant  gas  phase 

in  in  out  out  change 


Eq.  5.1.20 


where  the  energy  released  during  phase  change  is  a  function  of  the  mass  rate  of  fluid 
condensing  in  the  heat  exchanger  and  the  latent  heat  of  water: 

Eq.  5.1.21 

change 

Further  expansion  of  Eq.  5.1.20  yields: 


(mCT)  +(mCT) 

V  ^  ''main+ purge  V  ^  ^  cool 


-wLEq.  5.1.22 


where  the  main  flow  (combustion  products)  and  purge  flow  (air)  are  separated  in  the 
latter  half  of  the  equation.  The  main  flow  is  subdivided  into  nitrogen,  gaseous  water  and 
liquid  water  so  that  the  water  can  be  accounted  for  separately.  This  separation  is  possible 
because  the  experiment  was  always  run  with  an  equivalence  ratio  of  one  and  the  fuel  used 
was  hydrogen,  so  the  only  products  should  have  been  nitrogen  and  water. 

This  equation  is  manipulated  to  solve  for  Tavg,in  or  the  average  temperature  of  the 
gas  as  it  enters  the  heat  exchanger: 
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where 


Spurge  Q 
V  Ktal  ' 


Eq.  5.1.24 


It  should  be  noted  that  the  specific  heat  for  both  purge  and  main  flow  gases  vary 
significantly,  due  to  their  dependence  on  temperature.  For  this  reason  it  is  necessary  to 
estimate  the  purge  temperature  as  it  enters  the  heat  exchanger,  so  that  an  accurate  specific 
heat  may  be  calculated.  This  process  is  displayed  in  Appendix  C. 

The  mass  flow  rate  of  the  purge  and  main  flow  is  measured  in  the  purge  and  main 
manifolds  respectively.  It  is  assumed  that  the  mass  flows  are  conserved  and  do  not 
change  from  the  manifolds  to  the  exit  of  the  heat  exchanger.  The  gaseous  specific  heats 
are  attained  via  curve  fit  coefficients  (Turns  2006)  and  the  liquid  specific  heats  are  looked 
up  form  a  table  (Incropera,  et  al.  2007).  The  temperatures  for  the  exiting  gases  (purge  air, 
N2  and  H2O  (g))  are  measured  from  the  J-type  thermocouples  at  the  exit  of  the  heat 
exchanger  while  the  temperature  of  the  exiting  condensing  liquid  (H2O  (1))  is  assumed  to 
be  the  same  as  the  coolant  temperature  as  it  leaves  the  heat  exchanger. 

To  determine  the  mass  rate  of  liquid  water  condensing  in  the  heat  exchanger  an 
equilibrium  reaction  is  assumed  (Eq.  5.1.25). 

7/2  +^(6>2  +3.76A,)< - ^/72<^(g)  +  Eq.  5.1.25 
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This  assumption  lies  on  the  fact  that  the  reaction  is  taking  place  at  stoichiometric 
conditions,  which  is  believed  to  be  the  case.  It  is  also  assumed  that  the  purge  flow  is 
completely  composed  of  air  (O2+3.76N2).  Next  the  mole  fractions  for  the  main  exhaust 
products  are  determined  to  be: 

;f^,o=^  =  -347  Eq.  5.1.26 

Jv,=^  =  -653  Eq.  5.1.27 

Using  the  mole  fraction  of  water,  the  total  mass  flow  of  water  (gaseous  and  liquid)  that 
enters  the  heat  exchanger  may  be  calculated: 

^ H 2O (total)  ~  ^ main XH2O  Eq.  5.1.28 

Next,  steam  tables  (Sonntag  1991)  are  used  with  pressure  and  temperature  measurements 
of  the  gaseous  water  as  it  exits  the  heat  exchanger  to  approximate  its  density  (yO^^o(g))- 

This  density  is  used  to  determine  the  mass  flow  rate  of  gaseous  water  with  Eq.  5.1.29. 

m^pAV  Eq.  5.1.29 

The  area  is  based  off  of  the  inner  diameter  of  the  heat  exchanger  and  the  velocity  is 
determined  by  coherent  structure  velocimetry  provided  by  Schlieren.  Coherent  structure 
velocimetry  determines  the  velocity  of  the  large  and  small  scale  structures  present  in 
Schlieren  imagery.  It  accomplishes  this  by  tracking  the  structures,  depicted  by  the 
density  gradients,  from  one  frame  to  the  next.  The  frame  rate  and  the  frame  size  or  the 
spacial  resolution  (physical  area  represented  by  each  pixel)  are  both  known.  By 
calculating  the  distance  that  the  structure  travels  over  each  frame  and  dividing  by  the 
amount  of  time  elapsed  from  one  frame  to  the  next,  velocity  is  produced. 
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The  result  of  Eq.  5.1.29  is  the  mass  flow  rate  of  the  gaseous  H2O,  whieh  is 
subtraeted  from  the  total  mass  flow  rate  to  determine  the  rate  of  liquid  water  being 
produeed  at  the  exit  of  the  heat  exehanger: 

~  ^H^Oitotal)  ~^H20(g)  ^9-  5.1.30 

The  radiation  (qrad)  and  natural  conveetion  (qnc)  terms,  if  found  to  be  signifieant 
are  included.  These  calculations  can  be  found  in  Appendix  B.  With  this  information  it  is 
possible  to  return  to  Eq.  5.1.23  and  calculate  the  average  temperature  at  the  inlet  of  the 
heat  exchanger  (Tavg,in)-  It  is  necessary  to  iterate  on  this  temperature  because  the  inlet 
temperature  is  dependent  on  the  specific  heat  of  the  inlet  gases,  which  is  also  temperature 
dependent. 

When  the  analysis  was  performed  for  the  configuration  without  the  turbocharger  in 
the  flow,  the  result  of  Eq.  5.1.23  or  Eq.  4.2.3  was  considered  the  average  turbine  inlet 
temperature  and  when  performed  for  the  configuration  with  the  turbocharger,  turbine  exit 
temperature  was  calculated. 

The  inlet  temperature  calculated  during  these  experiments  using  Eq.  5.1.22  was  1139 
K  (1591  F)  with  the  turbocharger  and  1262  K  (1812  F)  without  the  turbocharger.  For  20 
Hz  at  a  fill  and  purge  fraction  of  0.5  and  0.5  respectively.  Looking  at  specific  work,  this 
translates  to  177  kW/kg  (1 .79  hp/lb/min).  It  is  possible  to  convert  the  specific  work  to 
average  work  by  multiplying  the  specific  work  by  the  mass  flow  rate.  Doing  so  yields 
4.45  kJ/s  (5.95  hp).  The  work  provided  in  the  initial  results  is  higher  than  expected.  This 
approach  may  provide  inflated  results  because  it  assumes  that  all  of  the  energy  removed 
from  the  flow  by  the  turbine  is  completely  converted  into  work.  In  reality  a  portion  of  the 
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energy  in  the  flow  is  heating  the  turbine.  The  turbine  housing  temperature  rose  above 
800  K  (1000  F)  during  extended  runs. 
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Appendix  F-  Uncertainty  Calculations 


The  first  sensor  to  be  examined  is  the  water  flow  meter.  The  water  flow  meter 
was  used  to  determine  the  mass  flow  rate  of  eoolant  through  the  heat  exchanger.  The 
meter  used  was  a  FTB4707  low  flow  liquid  flow  meter  was  used  which  has  an  accuracy 
of  ±  1.0%  (Omega,  Omega.com  2010).  The  temperature  of  the  coolant  was  measured 
using  T-type  Omega,  mini  quick  disconnect  thermocouple  probes,  both  at  the  entrance 
and  exit  of  the  heat  exchanger.  These  thermocouples  have  a  tolerance  of  1  K  or  0.75% 
(Omega,  ANSI  and  lEC  Color  Codes  for  Thermocouples,  Wire  and  Connectors  2010). 

Looking  at  the  heat  removed  from  the  system  by  the  coolant: 

^water  ~  ~  W  (^C^  ^ p,out'~^ out, coolant  )  ~  ^  Eq.  5.1.31 

Based  off  of  representative  values  a  nominal  heat  removed  by  water  is  calculated: 

q  =  (.347)[(4178)(310)-(4184)(292)] 

=  25487.8  W 


It  is  necessary  to  determine  all  of  the  terms  in  the  equation  before  solving: 

T  =  C,  r„„  -C,?;,  =[(4178)(310)-(4184)(292)]  =  73452 


5q 

dm 

dq 


out 


=  m 


C  -C 

y  out  in  J 


=  .347(4178-4184)  =  2.082 


.  =(.347)(.001)  =  ±.000347  kg/s 


w. 


Wj=±lK 


Now  applying  Eq.  4.2. 16  to  the  equation  for  heat  removed  by  the  water  yields: 

Wq  =  [(73452)"  (.000347)"  +(2.082)"  (l)" 

=  ±25.5  W  or  0.10% 
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The  air  mass  flow  rate  through  the  manifolds  is  governed  by  the  equation: 


m  =  c^Pq 


-  .0267  kg/s 


Eq.  5.1.32 


where  Pq  and  To  are  the  stagnation  pressure  and  temperature,  Cd  is  the  nozzle  diseharge 
eoeffieient,  D  is  the  nozzle  throat  diameter  of  the  converging-diverging  nozzle,  yis  the 
ratio  of  specific  heats,  R  is  the  gas  constant,  and  gc  is  the  gravitational  proportionality 
constant.  The  nozzle  throat  diameter  is  measured  to  the  nearest  .001  inches,  resulting  in 
an  uncertainty  of  ±.0005  inches  or  1.27  E'^  meters  (Engineering  2011).  The  nozzle 
throat  diameter  used  for  this  analysis  is  .252  inches  or  .0064  meters.  The  discharge 
coefficient  for  the  nozzle  is  0.991  with  an  uncertainty  of  ±1%.  The  thermocouples  used 
to  measure  the  temperature  in  both  manifolds  were  Omega  J-type  thermocouples  with  a 
published  accuracy  of  1  K  or  +  0.75%  (Omega,  ANSI  and  lEC  Color  Codes  for 
Thermocouples,  Wire  and  Connectors  2010).  Sensotec  model  TJE  pressure  transducers 
measure  the  pressure  in  each  of  the  manifolds  to  within  +  0.1%  accuracy  (Honeywell 
2005). 

Using  the  same  methodology  as  before,  the  components  of  the  equation  are: 
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1.4-1  j - ^ ^ - 

2(14-1)  I  (1.4)  (9 .8) 

2  J  '\|  (286,9)  (292) 


=  ,0273 


(,991)(.01)  =  -.00991 


Wp^  =  ±MPa 


Wp^  =  ±84  Pa 


-4.64  Wt^  =  ±1K 


(.0273)^  (.0099l)%(7.51  E"®)^  (84)^ 
-(8 .47)^  (1.27  £"')^+(-4.64  £’)^(-l)^ 


=  =.000117  kg/s  or  0.43 3% 


Next  the  purge  gas  temperature  as  it  enters  the  heat  exchanger  will  be  analyzed 
starting  with  the  Reynolds  Number.  The  tube  diameter  is  measured  to  the  nearest  1/16 
inches,  so  the  uncertainty  is  ±1/32  inches  or  0.00079  m: 
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Re  =  ^^  =  15294 
7iD|a 

5  Re  4 


dm  7iD|^  ;r(.055)(2.40xl0“^) 

5Re  -Am  (-4)(.0158) 


=  965467 


■  =  1111X1 


aD  ttDV  ;r (.055)  (2.40x10"') 

=  ±.0001 17  kg/s 

Wd  =  (.055) (.00079)  =  ±.000043  m 

Wr,  =  [(965467)"  (.0001 17)"  +(277212)"  (.000043)" 
=  ±113.587  or  0.74% 


ni/2 


Frietion  factor  uncertainty  is: 


/  =  (0.790Z«(Re)-1.64)'  =.028 
df  _  -3.20461 

5  Re  Re(L«  (Re) -2.07595) 
=±113.587 

(4.85  (113.587) 


=  -4.85  E 


-1 


Wf  = 


ni/2 


=  ±5.5  E-^  orO.lOVo 


Nusselt  number  (note:  this  formula  is  valid  for  a  wide  range  of  internal  flows)  uncertainty 


is: 


Nu  =  - 


(//8)(Re-1000)Pr 


=  44 


l  +  12.7(//8)‘'"(Pr"''-l) 

(.0139/‘^"  (Pr"''- 1)  +  .445)(Re- 1000)  Pr 


df 

dNu 


(r"(Pr- 

_  .0278/ Pr 

aRe~/''2(p//3_i)  +  .223 

=±5.5  E-^ 

=±113.587 


-l)  +  .223)" 
=  .002893 


=  1622 


^Nu  = 


(1622)"  (5.5  E-^  f  +  (.002893)"  (l  13.587)" 
=  ±.089  or  0.20% 


ni/2 
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Next  convection  coefficient  uncertainty  is  calculated: 


A  =  ^  =  27.6 
D 

dh  k  .037 

- =  —  = - =  .673  w^„  =  ±.089 

dNu  D  .055 

^ _ (44)(.037)  ^ =  ±0.00079  m 
dD  D"  .055" 


:  (.673)"  (.089)"  +  (538)"  (.00079)" 


=  ±.429  — —  or  1.56% 
m  K 


The  last  portion  of  the  purge  temperature  uncertainty  examines  the  purge  temperature 


formula  itself,  Eq.  5.1.19: 


T  =T  "NT  .-T  )  =  567.8  K 

purge  surf  \  surf  man  / 


Jlil  =  J  _1  0.305  rf  =  ±7.5  K 

^snrf 


f)T  ~\  — ^ 

^£^  =  ^^”"^^=0.701  w,_=±7.5K 


fTsu^-T...){PL) 


=  8.39  w,  =  ±4292 


=  -14681  =±.000117  kg/s 


(.305)"  (7.5)"  +(.701)"  (7.5)"  +(8.39)"  (.4292)"  +(14681)"  (.000117)" 


=  ±7.0  K  or  1.2% 


Moving  on  to  the  radiation  and  free  convection  of  the  heat  exchanger,  because  (3  for 


an  ideal  gas  is  only  a  function  of  temperature  (1 /temperature): 
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j3  =  —  =  .0033  K 
T 

avg 


-1 


1 


dj3 


w,,,^=±7.5K 


=.000011 


(.00001  if  (1)" 

=  ±.00008  K"'  or  2.5% 


1/2 


Ra 


gj0(T^-T^)D^  ^ 


37 


oa 


dRa  3gyg(T,-TjD^  3(9.8)(0.003)(324-292)0.168^  ^ 

5D  ua  (l.59  £-^)(2.25  E~^) 


Wj)  =  ±1.27  E  ^  m 


dRa  g/?D^ 


oa 

dRa  _  -gj3D^ 


oa 

dRa  _  g(T3-T^)D^ 
d/3  ~ 


=  389669 

=  -389669 


Wt,  =  ±7.5  K 


w,„=±lK 


=  4.16  E^ 


oa 


=±.00001  K 


-1 


nl/2 


(2.23  £^)"(l.27  E~^)\{2%9669)\l .5)  + 

(-389669)'  (-1)'  +(4.29  (.00001) 

=  ±2.94  .EV21.5% 

Next  the  uncertainty  of  the  Nusselt  number  for  a  long  horizontal  cylinder  is 
addressed: 
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2 


Nu=i 


dNu 

dRa 


0,60  + - 


1/6 


1  + (0.559 /Pr) 
0.061 


9/16 


8/27 


=  31 


=  6.01  £- 


m 


>v,,=±2.94  r 


= 


(6.01  (2.94 


1/2 


=  ±1.766  or  5.7% 


Once  the  uncertainty  of  the  Nusselt  number  is  known,  the  convection  coefficient 

uncertainty  may  be  determined: 

,  Nuk  ,  W 
h= - =  4.84 


D  mK 

dh  k  0.0263 


=  0.157  w^„=±1.76 


dNu  D  0.168 

dh  Nuk  5 

—  = - ^  =  28.9  Wn=±1.27fi  m 


5D 


= 


(0.157)"  (1.76)"  +(28.9)"  (l.27.E"^)" 


1/2 


=  ±.27 


W 
m  K 


or  5.7% 


The  uncertainty  for  the  radiation  and  free  convection  are  also  evaluated: 
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q„..=A;rDi(T.-T„)  =  503  (f 

%=!^  =  ;rDL(T.-T„)  =  104  >v,=±.27 

on 

dq 


m  K 


aD 

^^com 

dL 

^^com 

ai, 

^^com 

ar 


=  M(T^-T^)  =  2989  Wo  =±1.27  £-5  m 
=  hTvD  (T^  -T„  )  =  83  w^  =  ±.00079  w 
=  hTTDL  =  16  Wj^=±1.5K 


=  -hKDL  =  -16  Wt,  =  ±7.5  K 


w  = 

ofonv 


The 


(104)" (.27)"  +(2989)" (l.27  +(83)" (.00079)" 

+  (l6)"(7.5)"+(-16)"(-7.5)" 

=  ±172  IF  or  34% 

emissivity  of  6061-T6  aluminum  was  found  to  be  .02  with  an  uneertainty  of 

__  _  1  _ _ o  ^  N 


^rad  = 

=  13.7  W 


2 

II 

TiDLcriT^ -T^ 

M  =  684 

w  =  +.01 

ds 

\  s  00 

/ 

6 

^^rad  _ 

STrLcriT^-Tj 

)  =  81.4 

Wn  =  ±1.27  E~^  m 

dD 

\  s  00 

2 

II 

^;rD  cr(r/-F 

M  =  2.24 

Wj  =  ±.00079  m 

dL 

\  5  00 

/ 

^^rad  _ 

4f;7'DZ(Tr2  =  160  w.. 

=  ±7.5  K 

dT, 

s 

IS 

*l,.j  _ 

-AsttDLctTJ  = 

-106 

w.  =  ±7.5  K 

00 

Is 

(684)"  (.01)"  + 

(81.4)’ (1.: 

27  £^^)+(2.24)"(.C 

^qrad  = 

/  ^  ^  7  /  X  7 

/  ^  ^7 

/  X  7 

|l/2 


L+(160)  (7.5)  ■ 
=  ±1443  IF  or  10500% 
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Note  that  the  uncertainty  values  calculated  for  heat  transfer  only  reflect  the  error 
from  this  analysis.  Actual  uncertainty  values  for  heat  transfer  may  vary  over  a  much 
wider  range. 

The  formula  for  average  specific  heat  at  the  outlet  of  the  heat  exchanger  must  be 
analyzed  before  the  uncertainty  for  the  enthalpy  at  the  same  location  can  be  determined: 


fh  m  T 

^  _  purge  ^  ''‘'main  ^  _  i  i  qq  ^ 

^p,main+ purge  .  ^Ppurge  •  ^  Pmain  ^ 


total 


ac„ 


^Kain  Kta! 

ac„  1 

-c„ 


^  C„  =—!— (1320)  =  48889  w.  =±.Q00nikgls 

V  /  ^main  ^ 


.027 
1 


dm 


ffl 

purge  total 


mi 


(1007)  =  37296  _  =±.000117  kg !  s 


dCp  ^  Spurge  Q  ^main 

2  P  purge 


dm. 


ffl 

mam  total 


m 


C  = -42019 


2 


total 


mi 


mi 


=  ±.000234  kg ! s 

(48889)"  (.000117)"  +(37296)"  (.0001 17)"  +(-42019)"  (.000234)" 


w, 


c 

^  p  ,main+ purge 

out 


ni/2 


=  ±12.2 


kgK 


or  1.07% 


Next,  the  enthalpy  of  the  exhaust  gases  is  determined: 


r 


\ 

m  C  T 

total  p,avg  avg 

out  out  y 


=  (.027)(1139)(314)  =  9656  W 


^gas 
out 

^  -  8-48  Wq,  -  ±12.2 

^  =30.8  h.,=±2.2X 

01  out 


=  ±.000234  kg  / s 
J 


w, 


=  [(357646)"  (.000234)"  +(8.48)"  (12.2)"  +(30.8)"  (2.2)" 
=  ±150  W  or  1.5% 
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Knowing  the  uncertainty  of  the  mass  flow  for  both  the  main  and  purge  manifolds,  it  is 
possible  to  calculate  the  uncertainty  of  the  enthalpy  of  the  gas  as  it  exits  the  heat 
exchanger.  Again  J-type  thermocouples  are  used  for  the  temperature  measurement  of  the 
gas  at  the  exit  of  the  heat  exchanger. 

Now  that  all  of  the  input  arguments’  uncertainties  have  been  solved  for,  the  inlet 
temperature  uncertainty  can  be  considered: 


(T) 

main+purge 

in 

S(T)  main+purge 


\«.+qrad+q/c  +q  water 
out 


out  _ 


1 


1335  K 


main+purge 


dh 


gas 

out 


main+purge 


(.027)(1142) 


=  .032  w,  =±150  if 


3(T). 


main+purge 


S(T) 

main+purge 

in 

d(T) 

main+purge 


in 

1 


■  =  .032  w,  =±1443  1+ 

4ra(/ 


main+purge 


in 

1 


■  =  .032  w,  =±172  1+ 

Q/c 


main+purge 


dq 

a(T). 


main+purge 


main+purge 
in 

^gas  +^rad+<lfc+<l 


=  .032  =±25.5  W 

Hwater 


out  _ 


dm 


J  main+purge 


d{T)„ 


main+purge 


in 

hgas+(lrad+<lfc+<l. 


-49965  =±.000117  ytg/.s 


ac„ 


[^^p) 


^  =  1.02  Wc  =±12.2 


J 


main+purge 


kgK 


w,  = 


(0.032)"  (150)"  +(0.032)"  (1443)"  +  (0.032)"  (172)" 

+  (0.032)"  (25.5)"  +(-49965)"  (.000234)"  +(l.02)"  (12.2)" 
=  ±49  K  or  3.7% 


1/2 
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The  uncertainty  of  the  compressor  work  will  be  addressed  next: 


w  =  mC 

comp  p^ 


^  p  A  A  f  / 

Ti  ^  -Ti  =(.067)(1019)  (304) 

v^iy  j  V  V 


122803 

99739 


-(304)  =1270  W 


=±12.2  — — — 

Povg  ]^g 

w^=±\  K 


=  c  T.  ^  -T,  =71663  w.  =±.000117  kg/ 5 

aw  p 

um  y  J 


r  “il/2 

=  [(4.7)"  (12.2)"  +  (15)"  (1)"  +  (71663)"  (.0001 17)" 

=  60  WorA.1% 


Next  the  uncertainty  of  the  enthalpy  at  the  turbine  exit  will  be  addressed: 
^4  =wCpj4  =(.055)(1242)(1036)  =  70550  W 

^h^  =  jhT.=51  =±12.2^^ 

aCp^  '  kg  K 

dw 

— L  =  jnC=6S  Wj=±lK 

ax^  ■ 

^  =  C,T.  =1.28  E”  w,  =±.000117  kg/s 

dm 

r  -|l/2 

=  (57)"(l2.2)"+(68)"(l)"+(l.28  .E®)" (.0001 17)" 

=  714  r  or  1.0% 


The  uncertainty  of  the  corrected  turbine  inlet  temperature  was  calculated  as: 
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T  = 

3,„ 


Ko^p+Kud+K  ^  1270  + 1090  +  63990 
mC 

Pa 


(.055)(1233) 


=  978  K 


'{^comp+Kad+K) 


dC 


p* 


mC 


=  .79  =±12.2 


J 


Pavg 


Pa 


kgK 


ax. 


dW  mC 

comp  Pa 

ax. 


=  .015  w...  =±60  r 


^Kad 

ax. 


=  .015  w,  =±192  W 


dh, 


=  .015  w,  =±714r 


ax. 


'{^comp  ^  Kad  +  K) 


dm 


•  2/^ 


=  -17789  =±.000117 


Pi 


w  = 

■^corr 


ni/2 


(.79f  (I2.2f  +  (.015)"  (60)"  +(.015)"  (714)" 

+  (.015)"  (192)"  +(-17789)"  (.000117)" 

=  15  Kor\.5% 

Knowing  the  uncertainty  of  the  inlet  temperature,  the  work  performed  by  the  turbine 
may  be  assessed: 


-C,T,)  =  (.055)[(1302)(977)-(1294)(950)]  =  2351  W 


Sir, 

ax, 

corr 

ax. 


=  mC=12  w.  =±15  K 


Pz 


=  -mC„  =1\  =±4.7  K 

Pa  ^4 


dm 


W,i 


^  =  (CT.  -CX.)  =  42754  w.  =±-000117^ 

\  Pl  Pi  ^  J  *4 

=  [(72)"  (19)"  +(71)"  (4.7)"  +(42754)"  (.000117)" 


=  1130  IV or  48% 


Xhe  isentropic  relation  used  in  Eq.  4.2.13  is  probed  next: 
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r-i 


7;^  =  7;  (Pr)  r  =  (977)(2.4)i.4  =  1254  K 


dZ 


4s 


dZ 


Wt  = 


7~)  A 


v^3y 


r-i 


=  1.28  =±\S  K 


(1.28^(15)^ 


1/2 


=  19.2  K  or  5% 


Finally  the  efficiency  for  the  turbine: 


W.. 


Vt=- 


turb 


2356 


m(C^,T,_ ■  .055[(1302){977)-(1292)(951)] 

H  _  1 


=  .987 


0’^.*  (c„r3-c,„r„) 

brj, 


C  T 


dZ 


=  .000023  w,,  =±1130  IF 


=  .012  w.  =±15  7: 


It  _ 


T  W 

^4/^  turb 


As  ^Pi'l'icorr) 


=  .021  w.  =±19  7: 


dn  -W 

_6  ^ - - turb -  =±.0001177  kg/s 

dm  m^(C,3^3-C,4A) 


r  -|l/2 

=  (.000023)'  (1 130)'  +(.012)'  (15)'  +(.022)'  (19)'  +(18)'  (.0001 17)' 


=  .44  or  44% 
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